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Summary. — Simulating magnetic effects with cold gases of neutral atoms is a 
challenge. Since these atoms have no charge, one needs to create artificial gauge 
fields by taking advantage of the geometric phases that can result for instance from 
atom-light interaction. We review here some schemes that lead to the desired Hamil¬ 
tonians, either in a bulk geometry or in a lattice configuration. We also detail the 
relations between some general concepts of magnetism, such as gauge invariance, 
Landau levels, topological bands, and the features that can be generated in cold 
atoms setups. 
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Magnetic effects play an essential role in quantum physics. Notions or phenomena 
as diverse as gauge invariance, quantum Hall and Aharonov-Bohm effects, topological 
insulators, find their origin in the coupling between electromagnetic fields and particles 
carrying an electric charge. For a particle with mass M, charge q and velocity v, this 
interaction is described in terms of the Lorentz force F = qv x B or equivalently the 
Hamiltonian 


(i) 


H = 


(p - qMr)f_ 
2 M 
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where A is the vector potential associated to the field B and p the momentum operator 
of the particle. The study of magnetism constitutes an important aspect of the general 
quantum simulation program based on cold atomic gases, in which one hopes to emulate 
with these gases a large variety of quantum phenomena that one encounters in other 
fields of many-body physics, for example condensed matter physics [1]. However the fact 
that the atoms are neutral (q = 0) forces one to look for “substitutes” to real magnetism, 
for example using light beams with well chosen frequencies and directions. 

This course will constitute an introduction to the general subject of artificial (or 
synthetic) magnetic fields. Let us emphasize in advance that this subject is developing 
extremely fast and it is not possible in the limited space of these lecture notes to give an 
exhaustive presentation. Our goal is to present an overview of the domain, and orient 
the reader to more complete reviews such as [2, 3, 4] in order to deepen the subject. 

The outline of this lecture is the following. In the first section, we present some 
important features of magnetism such as gauge invariance and Landau spectrum. We 
discuss the Aharanov-Bohm effect, which will constitute later a guideline for identifying 
magnetic-like phenomena. We also present a first option to simulate magnetism, which 
consists in setting a gas of neutral particles in rotation. Sect. 2 is devoted to the sim¬ 
ulation of magnetism using geometrical phases, such as Berry’s phase. We show how 
to transpose this notion to particles moving in a light field and following adiabatically 
a dressed state, z.e., an eigenstate of the atom-light coupling. In sect. 3 we generalize 
the notion of geometric phase to the case of non-Abelian potentials, and discuss the 
particular case of spin-orbit coupling. Then sect. 4 is devoted to the important case of a 
particle moving in a lattice in the presence of a magnetic flux. The discretization of space 
introduced by the lattice considerably enriches the problem, with a single-particle energy 
spectrum that acquires a fractal nature. We then discuss various approaches to simu¬ 
late this lattice problem with neutral atoms, based either on optical lattice shaking or 
modulation (sect. 5) or on the use of a resonant coupling between atomic internal states 
(sect. 6). Finally we conclude by outlining some perspectives of this field of research, 
when one tries to go beyond the single-particle problem considered in these notes. The 
first appendix attached to this chapter discusses the eigenvalue problem for a particle 
in a uniform magnetic field (Landau levels). The second appendix details the notion of 
topology for a particle in a lattice, and outlines the characterization of these topological 
properties in terms of an integer number, the Chern index. 

1. — Magnetism and quantum physics 

IT. Gauge invariance. - We will essentially consider in these notes magnetic effects 
induced by a static magnetic field B{r). This magnetic field is divergence-free, = 0, 
which expresses the fact that the flux of B across any closed surface is zero or, in other 
words, that there does not exist free magnetic monopoles. From this constraint, one 
deduces that B can be written as B = V x A, where A is a vector field. It is clear that 
this writing is not unique since two vector potentials A and A' such that V x {A-A) = 0 
correspond to the same magnetic field and have the same physical effects. More precisely, 
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since a vector field with zero curl can be written as the gradient of a scalar function y(r), 
the set of vector potentials that correspond to a given magnetic field B(r) are related to 
each other by 

(2) A\r) = A(r)+ Vx(r), 


where x( r ) is an arbitrary (sufficiently regular) function. 

Consider for example the case of a uniform field B = B u z , where Uj denotes the unit 
vector along the j axis. Among an infinite number of options, three convenient choices 
for A are the following vector fields in the xy plane: 

(3) A(r) = Bxu y , A(r) = —Byu x , A(r) = — (xu y - yu x ), 

where the two first expressions correspond to the Landau gauge and the third one to the 
symmetric gauge. The function y allowing one to change from one gauge to the other is 
in this case proportional to Bxy. We note that these gauge choices break the in-plane 
translational and rotational invariances of the initial problem. 

We now turn to the description of the motion of a charged particle in the magnetic 
field B within Lagrangian and Hamiltonian formalisms. Imposing the equation of motion 
mr = qr x B , one finds that the class of suitable Lagrange functions is 

(4) L(r, f) = ^Mr 2 + q r • A{r). 

The Euler-Lagrange equations 


( 5 ) 


dL_ _ _d fdL\ 
dn d t\dri)' 


ri = x,y,z, 


then provide the desired result. The Lagrange function (4) clearly depends on the gauge 
choice. However one can check from eq. (2) that Lagrange functions corresponding to 
gauge choices for the same B field differ only by a total derivative with respect to time, 
hence do correspond to the same physical problem. 

The Hamilton function H(r,p) associated to the Lagrange function L(r, r) is obtained 
by calculating first the conjugate momentum 


(6) p = Vf.L(r, r) = Mir + qA(r ), 

and then making the Legendre transform H(r,p) = p r — L(r,f), which leads to the 
result (1). When turning to quantum mechanics, canonical quantization amounts to 
associate operators r and p to the position and canonical momentum of the particle, 
with the quantization rule [rj,pk] = i hSj^ [5, 6]. We will always take in the following 
the standard choice p = —i^V r , which ensures that these commutation relations are 
fulfilled. 
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It is clear that if a given wavefunction ^(r, £) is a solution of the Schrodinger equation 

(7) = ( - mV 2 -^ (, ' )) % (r, t ), 

it will generally not be a solution of the Schrodinger equation for another gauge choice 
A'(r) that is deduced from A(r) by the gauge change of eq. (2). The remedy to this 
problem is simple: one has to impose that the wavefunction is also modified in a gauge 
change. More precisely a gauge transformation in quantum mechanics is defined as the 
simultaneous substitutions 

(8) A(r) —» A'(r) = A(r) + Vy(r), 

(9) ^(r,^)—> — 7V(r,t), with T = exp[igx(r)//i]. 


One can then check that if ip is a solution of the Schrodinger equation (7) for the vector 
potential A, then if' is a solution of the Schrodinger equation for A'. 

1'2. Cyclotron motion and Landau levels. - Consider a charged particle placed in a 
uniform magnetic field B parallel to the £ axis. We will restrict in the following to 
the dynamics of the particle in the xy plane. If the motion of the particle is described 
by classical mechanics, the particle undergoes a uniform, circular motion with angular 
frequency 


( 10 ) 



called the cyclotron frequency. In quantum mechanics, dimensional analysis indicates 
that a natural energy scale appears, huj c , as well as the length scale 



called the magnetic length. For an electron in a field of 1 T, the cyclotron frequency and 
the magnetic length take the value uj c /2tt = 28 GHz and ^ mag = 26 nm. 

The reason for which the magnetic length scale appears physically in the quantum 
problem can be understood as a consequence of Heisenberg inequality. Classically, the 
size ro and the velocity vo of a cyclotron orbit are linked by the linear relation vo = ^cTo- 
the smaller the orbit ro, the smaller the corresponding velocity vq. In quantum terms 
one cannot prepare the particle in a state where both the position and the velocity are 
arbitrarily well known and Ar Av > h/(2M). The magnetic length can be understood 
as the minimal cyclotron orbit size compatible with this inequality. 

The energy spectrum of a charged particle in a uniform B field is remarkably simple 
and consists of equidistant Landau levels: 



( 12 ) 


n G N. 
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To prove this result we introduce the kinetic momentum 
(13) II = p — qA(r ) 

and rewrite the Hamiltonian (1) as 

< 14 > Jj = ^(ni + n;). 

Contrarily to r and p, two components of II do not commute but their commutator is a 
constant for a uniform field: 

(15) [tl x ,ily] = ihqB. 

We are then facing a problem that is formally equivalent to the search of the eigenvalues 
of a harmonic oscillator Hamiltonian (P 2 + X 2 )/2 with [X,P\ = i, hence the structure 
in equidistant energy levels. 

The reasoning above is gauge-independent; to go further and determine a basis of 
eigenstates of the Hamiltonian, one needs to specify a gauge choice. The procedure is 
outlined in Appendix 1 both for the Landau gauge and the symmetric gauge. One finds 
that each Landau level has a macroscopic degeneracy Af , proportional to the area A 
accessible to the particle in the xy plane, A f = A/2 / irf% riaig . This result can be interpreted 
as the fact that each independent state in a given Landau level occupies the area 27r^ ag . 
Since ^ ag l/^> this degeneracy can also be written A f = where we have 

introduced the flux = AB of the magnetic field through the accessible area and the 
flux quantum <f>o = h/q. 

1'3. The Aharonov-Bohm effect. - In their famous 1959 paper [7], Aharonov and 
Bohm( 1 ) proposed a gedanken experiment that illustrates a remarkable feature of quan¬ 
tum mechanics: One can detect the presence of a magnetic field using measurements 
made on particles that have never penetrated the regions of non-zero field. 

The geometry proposed by Aharonov and Bohm is illustrated in fig. 1. A two-path 
interferometer is pierced by an infinite solenoid, which is surrounded by a potential 
barrier so that the particles cannot enter inside. When a current runs in the solenoid, a 
non-zero magnetic field appears inside the solenoid, but the field is always zero outside. 
However the fact of running a current through the solenoid causes a displacement of the 
interference pattern detected on the screen. This displacement corresponds to a phase 
shift of one arm of the interferometer with respect to the other one 

(16) Atp = 2 tt L, 

u>0 


( 1 ) We use here the standard terminology for this effect, although a very similar discussion was 
made 10 years before by Ehrenberg and Siday [8]. 
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Fig. 1. - Gedanken experiment of Aharonov and Bohm [7]. A non-zero magnetic field inside the 
solenoid causes a displacement of the interference pattern detected on the screen. 


where <f> is the magnetic flux through the solenoid and <f>o = h/q the flux quantum. 

The Aharonov-Bohm phase is a geometric phase, in the sense that it does not depend 
on the velocity of the particle, nor on the time it takes to go from the emitter to the 
detection screen. It can even be called a topological phase in this particular geometry, 
since one can deform at will the trajectories without changing A<p, as long as one keeps 
one path on each side of the solenoid. 

The phase factor e lA( ^ is gauge invariant since it is expressed in terms of the magnetic 
field itself [9]. If one insists in expressing Acp in terms of ‘local’ quantities, z.e., quantities 
defined in the region accessible to the particle, then it has to be written in terms of the 
vector potential: 

(17) Aif = JJ B z (x, y) dr Ay = f jf qA(r ) • dr, 

where the line integral is taken along the oriented contour C consisting of the two paths 
of the interferometer. 

We briefly outline below one possible reasoning leading to eqs. (16) and (17). The 
idea is to choose different gauges on different patches of space, and then connect these 
patches together [7]. In the absence of current, we can make the gauge choice A = 0 in 
the whole space. When a current circulates in the solenoid, A(r) is non-zero outside the 
solenoid but B remains null in this outer region, hence we have 

(18) V x A(r) = 0 out of the solenoid. 

Once a choice for A(r) has been made over the whole space, consider the two simply 
connected regions (z.e., regions without a hole) sketched in fig. 2, which we label I and 
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Fig. 2. - Aharonov-Bohm geometry seen from above. The light grey zones are two patches over 
which the vector potential in the presence of current can be written as A(r) = V%(r). The 
function \ is not the same for the two patches. 


II. Over each region we can unambiguously define the two scalar functions xi and Xu by 

(19) Xi.nM = [ A(r')-dr', 

Jo 

where the line integral goes from the emitter to a point in region I or II without leaving 
this region. Since patches I and II are simply connected, the functions Xi,n( r ) are indeed 
single-valued. 

Consider the matter wave ijj\ that propagates in the region on the left of the solenoid, 
inside patch I. For this wave, switching on the current in the solenoid corresponds to the 
modification 

(20) A(r) = 0 —» A(r) = Vxi(r), 

which can be viewed as a gauge change. This gauge change must be accompanied by the 
modification of the wavefunction (9). The same arguments holds for the matter wave 
propagating in the region on the right of the solenoid, inside patch II. Matter-wave 
interference at a point r on the detection screen is related to 

(21) tl>*(r)ip T (r) = exp([i«?(xn(r) -%i(0)A] d° } (OV’FW* 

where ib^j are the wavefunctions in the absence of current. The phase entering the 
prefactor can be written 

l>r nr 

(22) XII (r)-Xi(r)= A(r) ■ dr- A(r)-dr, 

J 0,Cn J 0,Cn 

where C\ and C\\ are two paths going from 0 to r and located respectively in patches I 
and II. This corresponds to the announced result (17). 

The geometrical aspect of the Aharonov-Bohm phase is a hallmark of magnetism in 
quantum mechanics [10]. The general program of magnetism simulation can be viewed 
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as a search for possible ways to induce a geometric phase with similar properties on a 
neutral particle. 

1*4. Rotating gases. - Among the various methods that have been developed in order 
to simulate magnetism, rotating the trap holding the particles is probably the simplest 
one from a conceptual point of view. The idea is to take advantage of the similarity 
between the magnetic Lorentz force, qv x £?, and the Coriolis force which appears in the 
frame rotating at angular velocity £4 around the z axis: 

(23) ^Coriolis = 2Mj) X Cl, Ct = ftu z . 

More precisely the Hamiltonian in the rotating frame is 

(24) H = ^ + V(r) - nL z 

where V (r) is the trapping potential in the xy plane and L z = xp y — yp x is the component 
of the angular momentum operator along the z direction. This Hamiltonian can be 
rewritten in the more suggestive way 

(25) H = + V(r) + Rentrif. (f) 

where the effective vector potential and the corresponding magnetic field are 

(26) qA(r) = M£l (xu y — yu x ), qB = 2Mf2, 

justifying the analogy between the Coriolis and the Lorentz forces. In addition to the 
Coriolis force, it is well known that the centrifugal force MQ 2 r also appears in the ro¬ 
tating frame. It is indeed present in the Hamiltonian (25), since the centrifugal potential 
is given by 

(27) VentrifV) = 

To implement this method in practice, one can use the harmonic trapping potential: 

(28) V(r) = V 0 (r) + 5V(r), V 0 (r) = ^Mlj 2 (x 2 + y 2 ), SV(r) = ^Mu 2 (x 2 - y 2 ). 

Here the coordinates x, y refer to basis vectors u Xl u y that rotate around 2 with frequency 
ft with respect to an inertial frame of reference, and the dimensionless parameter e 
characterizes the strength of the anisotropy. The success of this method can be tested 
by implementing it with a superfluid sample, and looking for the vortex lattice that is 
expected in steady state (see e.g. [11, 12, 13]). 



Introduction to the physics of artificial gauge fields 


9 


It is interesting to see how one can implement exactly the magnetic Hamiltonian 
(1) with this rotation method. Comparing eq. (1) with eq. (25), we see that two con¬ 
ditions need to be simultaneously fulfilled: (i) The centrifugal potential has to balance 
the trapping potential, so that Q = cv, which means that one has to resort to an extra 
potential (possibly quartic like in [14]) in order to confine the particles; (ii) The trapping 
anisotropy e must tend to zero. At first sight the latter condition seems impossible to 
fulfill, since nothing sets the gas in rotation in this case. However this difficulty can be 
circumvented by using the evaporative spin-up technique invented by the Boulder group 
[15]. It consists in preparing first a cloud rotating at a moderate angular frequency (for 
example Q ~ 0.7 u) using a non-zero e, switching off this anisotropy, and then performing 
a selective evaporation that eliminates atoms with an angular momentum below average. 
In this way the angular momentum per remaining atom increases in the course of evap¬ 
oration. At the end of the evaporation, the search for the equilibrium state of the gas 
must take into account the angular momentum L z of the trapped particles, which can be 
done using the Lagrange multiplier method. One then ends up by minimizing the same 
mathematical quantity as in eq. (24), where Q now stands for the Lagrange parameter 
associated to L z . Experimentally the Boulder group could reach Q = 0.993 uj with this 
method [15]. 

2. — Geometric phases and gauge fields for free atoms 

The notion of adiabatic evolution is frequently encountered in classical and quantum 
physics, when one looks for the evolution of a system whose state depends on external 
parameters. Let us consider this type of evolution and focus on the case where the 
external parameters travel along a closed trajectory. This means that these parameters, 
which we represent by a vector A, evolve slowly from time 0 to time T, and take at time 
T the same value as at initial time: 

(29) A(0) —> A {t) —> A (T) = A(0). 

It is then natural to assume that the system under study is in the same state at time T 
as at time 0. However this is not always true, as exemplified by the well-known Foucault 
pendulum. Suppose that the pendulum oscillates linearly at a initial time. Twenty-four 
hours later, the Earth has made a full rotation around its axis; hence the suspension 
point of the pendulum, playing here the role of the external parameter, is back to the 
same place. However the oscillation plane of the pendulum has rotated by an angle that 
depends on the latitude at which the pendulum is located. 

The goal of this section is to take advantage of this anholonomy in the quantum case, 
in order to generate geometric phases (Berry’s phase) that can mimic the Aharonov- 
Bohm phase even if the considered particle has no electric charge. 

2T. Berry’s phase. - Consider a quantum system whose Hamiltonian depends on a 
continuous parameter A. We assume for the moment that this parameter is controlled by 
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an external operator. The parameter A can stand for a real number or a set of numbers. 
For each value of A we denote by (^(A)) the eigenstates of the Hamiltonian and E n (X) 
the corresponding energies: 

(30) H( A)|V’n(A)) = E n (X) IV’n(A)). 

We suppose that the set (| , 0 n (A))} forms an orthonormal basis of the Hilbert space for 
each value of A. We are interested in the evolution of the state vector of the system 

(31) |W)) = ^2 c n(t) ivwvo]) 

n 

when A slowly evolves in time. More specifically we suppose that the system is prepared 
at initial time in one particular eigenstate {'ipg): 

(32) q( 0) = 1, c n (0) = 0 if n 7 ^, 

and we suppose that the evolution of A is slow enough for the adiabatic theorem to hold 
[16]. In first approximation the state of the system is thus at any time t proportional 
to |'0r[A(t)]), and the evolution of the coefficient cg{t) is obtained from the Schrodinger 
equation: 


(33) 


i h eg — 


Eg{t) - \h\ • (i/jelVi/je) 




Eg{t) - A • Ag(X) 


eg. 


Here we introduced the real vector called Berry’s connection : 


(34) Ag(X) = ih(^g\V^g), 

which will play the role of the vector potential (more precisely of qA) in the following. 

Suppose now that the parameter A follows a closed contour C in parameter space, so 
that A (T) = A(0). The integration of the evolution equation (33) is straightforward and 
gives: 

(35) c £ (T)®e i$dyn ( T )e i$ge ° m ( T ) Q ( 0) 


where we introduced the dynamical phase 


(36) $ d y n (T) 

and the geometrical phase [17] 


1 

h 


E((t) d t 


(37) 


<|)geom. _ y 

h 



X • Ag [A] d t 


1 

h 


Ag( A) -dA. 
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The dynamical phase is the usual phase that appears also for a time-independent problem. 
The geometrical phase (Berry’s phase) only depends on the “trajectory” of the parameter 
A during the evolution, and it is independent of the duration T. Both terms exp(i<F dyn -) 
and exp(i < h geom -) are physical quantities in the sense that they are gauge-invariant: they 
are unchanged if one modifies the definition of the eigenstates ^[A]) by multiplying 
them by an arbitrary phase factor. 

In the following we will restrict our discussion to the case where the parameter A 
evolves in a two-dimensional or a three-dimensional space. For example this parameter 
can stand for the position of a particle or for its quasi-momentum when it moves in a 
periodic potential. We can then introduce Berry’s curvature , which plays a role similar 
to a magnetic field: 

(38) Bt = V x A*. 

This quantity is a real, gauge-invariant, vector field. In full analogy with the Aharonov- 
Bohm phase we can rewrite the geometrical phase accumulated by the particle when the 
parameter A moves along the closed contour C : 

(39) $ge° m (T) = j i JJ Be - d 2 s 

where <S is a surface delimited by the contour C. 

2 *2. Adiabatic following of a dressed state. - The notion of geometric phase can be 
directly adapted to the case of an atom slowly moving in a monochromatic laser field [18]. 
Two types of degree-of-freedom come into play. First the center-of-mass motion of the 
atom can be described in terms of its position operator r and momentum operator p = 
—iftV r . Second the internal dynamics of the atom corresponds to transitions between 
electronic states induced by the laser field. Within the rotating-wave approximation, this 
internal dynamics is described by a time-independent Hamiltonian iA nt (r), so that the 
total Hamiltonian reads: 

(40) Htot = ^ + ffi„t(r). 

Treating first r as an external parameter, we define the dressed states as the eigen¬ 
states of the internal Hamiltonian 

(41) H int (r)\ip n (r)) = E n (r) |i(r)). 

At any point r, the set {| 'ip n { r )}} forms a basis set for the Hilbert space associated to 
the internal degrees of freedom of the atom. Then we consider the total quantum state 
of the atom and write it as: 


4 = ^2<l> n (r,t)\ip n (r)). 


(42) 


n 
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The characterization of the atom dynamics amounts to determining the probability am¬ 
plitudes 0 n (r, t) to find the atom at point r in the internal state \^ n {r)). We now assume 
that the particle is prepared at initial time in a given dressed state |^). We also suppose 
that it moves slowly enough so that it essentially remains in this state at any time and 
the contribution of the V’n’s for n y £ can be neglected. The validity of this assumption 
will be discussed in subsect. 2‘4. Starting from the Schrodinger equation 

dtir / h 2 A \ 

(43) ih— = H tot V(r,t) = + 


we can obtain an equation for the relevant probability amplitude fa. After some algebra, 
the result can be written in the form: 


(44) 


i h 


dfa 

dt 


(p-Aj(r)y 
2 M 


+ Ei(r) + Ve(r) 




It has exactly the structure of a scalar Schrodinger equation for a charged particle (with 
q = 1 by convention) moving in the magnetic field associated to the vector potential 

(45) Ag{r) = i h (fa\Vfa), 

which is nothing but Berry’s connection introduced in eq. (34). In addition the particle 
feels a potential that is the sum of two terms. The first one is simply the energy Eg(r) 
of the occupied dressed state, and the second one is the additional scalar potential 

(46) W-teEl< V V^>| 2 . 

n^i 

Physically Vg(r) represents the kinetic energy associated to the micromotion of the atom, 
as it makes virtual transitions between the effectively occupied dressed state V’r and all 
other dressed states n ^ £ [19, 20]. 

2‘3. The two-level case. - For simplicity we now focus on the case of an atom with only 
two relevant internal states. We neglect any process related to the spontaneous emission 
of photons. This situation can occur if we are dealing with the intercombination line 
of an atom with two outer electrons, like magnesium, strontium or ytterbium (fig. 3a). 
In this case the radiative lifetime of the excited atomic state is much longer than the 
relevant experimental time scales, and the atom-laser coupling is characterized only by 
the Rabi frequency k and the detuning A between the laser and the atomic frequencies. 

The two-level formalism can also be applied to the case of a A level scheme (fig. 3b), 
with a Raman transition between two sublevels of the atomic ground state. This scheme is 
relevant for alkali-metal species as well as erbium or dysprosium. Then if the detuning A e 
is large enough, the excited state e that plays the role of a relay for the Raman transition 
can be formally eliminated from the dynamics. One is then left with a two-level problem 
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Fig. 3. - (a) Quasi-resonant excitation of a two-level atom, (b) Stimulated Raman transition 
between two sublevels of the atomic ground state. 


in the subspace {|#i), |# 2 )}, with again two parameters: the Raman detuning A and the 
two-photon Rabi frequency k: 


(47) 


K = 




* 

b 




where Kj is the one-photon Rabi frequency for the transition gj 0 e. 

In both cases the atom-laser coupling, i.e., the Hamiltonian H[ nt responsible for the 
internal dynamics of the atom, can be written as a 2 x 2 matrix( 2 ) in the basis {| g), | e)} 
or {Iffi), \g 2 )}- 


(48) 


Vint = 


h 

2 




where both the Rabi frequency k and the detuning A may depend on the atom center- 
of-mass position. We define the generalized Rabi frequency £2, the mixing angle 0 and 
the phase angle </>: 


(49) = (A 2 + l^l 2 ) 1 / 2 , cos 0 = , sin#=^, k = \k\ , 

so that the atom-laser coupling becomes: 


(50) 


_ M2 / cos 6 
“ ~2~ Ve i0 sin(9 


e 1( ^ sin Q\ hQ 
-cos 0 ) = ~2 n a ' 


( 2 ) For the A scheme, the levels g± and #2 are also light-shifted by the lasers a and b. For the 
laser coupling {h/ 2) (K a \e){gi \ + ^b\e)(g 2 \) + H.c., these light-shifts are h\Kj | 2 /4A e , j = a,b. We 
assume that these light-shifts have been reincorporated into the definition of A. 
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where n is the unit vector characterized by the spherical angles ( 0 , 0 ) and the frj 
(j = x,y,z) are the Pauli matrices. The eigenvalues of this matrix are ±hQ /2 and 
the corresponding eigenstates read: 


(51) 


IV+) = 


/ cos( 0 / 2 ) \ 

Ve^sin(0/2)J ’ 


IV’-) 


/—e 1 ^ > sin( 6 ) / 2 )\ 

\ cos( 0 / 2 ) y * 


Suppose now that the state of the atom follows adiabatically one of the dressed states 
10±). Berry’s connection and curvature associated to this state are: 

(52) A± = (cos0 — 1 ) V0, 

(53) £± = ±^V(cos 6 >) x VA 
and the scalar potential reads: 

(54) v± (r ) = m [(' ve)2 + sin2 6 {y(t))2 . • 

Since we are interested here in the simulation of magnetism, we look for a non-zero 
Berry’s curvature B ±, which can be achieved only in the presence of non-zero gradients 
of the phase 0 as well as the mixing angle 0. The latter can be obtained either via a 
gradient of intensity (Vft) or a gradient of detuning (VA). 

Before going further we note that we can perform a gauge transformation in this prob¬ 
lem by multiplying the expression (51) of the dressed states by an arbitrary phase factor. 
In particular one can choose instead of eq. (51) the following dressed-state definition 


(55) 


W + ) = 


fe 1 ( ^cos((9/2)\ 

V sin(6>/2) J ’ 


W-) 


( — sm(0/2) \ 
cos(#/ 2 ) J ’ 


leading to Berry’s connection 

(56) A'± = ±^ (cos 0+1) VA 


Berry’s curvature 15± is unchanged in this procedure since it is a gauge-invariant quantity. 

As an example consider the case of the A scheme of fig. 3b, where the two laser beams 
are supposed to be plane waves propagating along the direction, with wave vectors 
k a « —kb = ku x . The two-photon Rabi frequency can thus be written K>(r) = e 2lb , 
with the constant amplitude k,q and the phase angle 0 =@ 2 kx. We also suppose that 
the detuning varies linearly in space along the y direction: A(r) = A'y, this can be 
achieved for instance using the differential Zeeman shift between the sublevels g a and g 5 
in a magnetic field gradient. A natural length scale appears in the problem 


( 57 ) 


£ = k 0 /A, 
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Fig. 4. - Grey scale representation of the magnitude of Berry’s connection generated in a A level 
scheme (see text). The Aharonov-Bohm-Berry phase for the contour C is ~ 27r. 


which represents the width of the region centered on the x axis where the Raman exci¬ 
tation is quasi-resonant (Rabi frequency larger than detuning). This geometry provides 
the necessary ingredients for a non-zero Berry curvature, z.e., non-zero and non-collinear 
gradients for the two angles </> and 0. It leads to 

(58) B±(r) = ±B 0 C 3 '\y) u z 
where 

hk 1 

(59) *0 = T , C { y) = TW - 

This artificial magnetic field is translationally invariant along the x axis. It is maximal 
for y = 0 and decreases as \y\~ 3 when y —> oo. Let us comment briefly on its magnitude 
Bq. Consider the rectangular contour of fig. 4, of extension \ = 2ir/k along x and i along 
t/, inside which the field is approximately uniform. The Aharonov-B ohm-Berry phase 
associated to this contour is 

(60) $ geom ' = \JJ & ' « d 2 r « 2 tt. 

This means that in this geometry, a superfluid gas should exhibits vortices essentially lo¬ 
calized around the line y = 0, with a typical linear density of one vortex every wavelength 
A along the x axis. 

A scheme very similar to what we just described was implemented at NIST in 2009 
in the group of I. Spielman [21]. The non-zero value of A' was achieved with a magnetic 
gradient of a few gauss per centimeters, and a value of £ of a few tens of optical wavelength 
A. When placing a rubidium Bose-Einstein condensate in this configuration, the NIST 
group observed the nucleation of quantized vortices, proving thus the existence of a 
non-zero artificial gauge field. 
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2*4. Validity of the adiabatic approximation. - In what precedes we have assumed that 
the two-level atom prepared in a given internal dressed state could follow adiabatically 
the dressed state when moving in the light field. Here we briefly discuss the validity of this 
approximation. Generally speaking the validity criterion for the adiabatic approximation 
is that the angular velocity of the eigenstate should be much smaller than all relevant 
Bohr frequencies involving this eigenstate (for more details, we refer the reader to [16]). 
Here, for an atom moving at velocity v, the angular velocity for a dressed state is ^ 
kv , since the expressions (51) of the dressed states vary typically on the length scale 
k~ l = A/27r. For a two-level atom, the relevant Bohr frequency is the generalized Rabi 
frequency Vt so that the validity criterion reads kv <C Vt. Since the relevant velocities 
for an atom absorbing and emitting photons from counter-propagating light waves are 
at least of the order of the recoil velocity hk/M , we find as a necessary condition 

(61) E r < m, 

where we introduced the recoil energy E r = h 2 k 2 /2M. In practice, eq. (61) constitutes 
a relevant criterion to check the applicability of the adiabatic Schrodinger-type equation 
(44). 

The recoil energy also enters when one evaluates the maximum Berry’s curvature that 
can be created with laser light. Assuming that both the phase and the mixing angles 
vary significantly over the reduced wavelength k~ l = A/27r, we find from eq. (53) that 
the corresponding cyclotron frequency uj c = B/M is such that 

(62) hhj c ~ E r . 


In order to reach such large values, a natural strategy is to switch to optical lattice 
configurations, in which case atom-laser couplings indeed vary significantly over A/27r. 
This will addressed in sect. 4. 


2*5. Spontaneous emission and recoil heating. - The use of a Raman transition between 
sublevels of the electronic ground state of an atom is an efficient way to reduce the 
heating originating from the random momentum recoils caused by spontaneous emission 
processes. However the heating rate may not always be lowered down to an acceptable 
level, especially for species from the alkali-metal series. Another more favorable class of 
atoms is the lanthanide family, with species like erbium or dysprosium that have recently 
been brought to quantum degeneracy (see e.g. [22, 23]). To facilitate the comparison 
between species, we define the dimensionless merit factor 


(63) 


^eff 

7 


where 7 is the spontaneous emission rate, and we analyze how it can be maximized. Our 
discussion follows closely that of refs. [24, 25]. 
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ns 


Fig. 5. - Left: Resonance transition ns CA np of an alkali-metal atom. The fine structure of the 
excited state leads to a splitting of the resonance line into two components D i and D 2 . Right: 
hyperfine structure of the ground state represented for a nuclear spin I = 3/2 (case of 7 Li, 23 Na, 
39 K, 41 K, 87 Rb). The splitting between the various Zeeman sublevels is produced by an external 
magnetic field. 


Alkali-metal species. Consider an alkali-metal atom with a electronic ground state | g) = 

|ns), irradiated by a monochromatic light beam with a frequency close to the resonance 
transition ns 0 np. Because of fine-structure coupling, the first excited electronic level is 
split into two sub-levels npi/ 2 and np 3 / 2 , giving rise to two resonance lines D\ and D 2 . In 
addition, because of the non-zero value of the spin I of the nucleus, the ground state ns 
and the excited states np 1 / 2 and np 3 / 2 are also split by hyperfine interaction. In particular 
the ground state is split into two sublevels with angular momentum F == I d= 1/2 (fig. 
5). The degeneracy between the various Zeeman states can then be lifted by an external 
magnetic field and resonant Raman transitions may occur between adjacent Zeeman 
states when one irradiates the atom with light beams with well chosen frequencies and 
polarizations (figure 5, right): | g,F,mp) ^ |g,F ,mp ± 1). Here we will assume for 
simplicity that the relevant detunings are large compared to the hyperfine splitting, and 
treat the atomic transition as if the nuclear magnetic moment were zero. 

Taking into account the Clebsch-Gordan coefficients for J g = 1/2 J e = 1/2 and 
J g = 1/2 J e = 3/2 transitions, the corresponding matrix element is (up to a multi¬ 
plicative factor of order unity that depends on the choice of polarization and hyperfine 
structure) 



Here Ai = uo — and A 2 = uo — uo 2 are the detunings of the laser with respect to the 
resonance transitions D\ et D 2 . The contributions of these two transitions appear with 
opposite signs in (64), because of a ‘destructive’ interference between the two paths going 
from \g,F, mp) to \g,F, mp d= 1), and passing via npi/ 2 or np 3 / 2 . This interference is 
linked to the electric-dipole character of the atom-light coupling, qr ■ £, where r is the 
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position of the outer electron and £ the electric field of the light. This coupling acts 
on the orbital degrees of freedom of the electron, but not on its spin. The fact that the 
electron spin can nevertheless be changed in a Raman transition originates from the spin- 
orbit coupling, which is responsible for the lift of degeneracy between npi/ 2 and np 3 / 2 . 
Now, when the detuning of the light is large compared to the fine structure splitting 
(A e = Ai ~ A 2 ), the effect of spin-orbit coupling becomes negligible and the Raman 
coupling amplitude rapidly vanishes (like 1/A 2 ), faster than what one could have naively 
expected from (47). 

The spontaneous emission rate is (again up to a multiplicative factor of order unity): 


(65) 


Tk 2 f 1 2 \ 

7 ~ — \3Af + 3A|J ’ 


where T _1 stands for the radiative lifetime of the excited states npi/ 2 and np 3 / 2 . The 
Raman matrix element and the spontaneous emission rate are thus both proportional to 
the light intensity (ft 2 ), so that the merit factor (63) does not depend on this intensity but 
only on detuning. The optimal merit factor is obtained by taking Ai « — A 2 « Af. s ./2, 
where A s .f. = Ai — A 2 is the fine structure splitting between npi/ 2 and np 3 / 2 : 

(66) M = |^. 

For rubidium atoms, the merit factor is M. = 8 x 10 5 . Choosing the effective Rabi 
frequency ft e ff equal to the recoil frequency cj r = E r /h = 27r x 3.8 kHz, we infer the 
photon scattering rate 7 = 0.03 s -1 . In average, the random recoil associated with the 
spontaneous emission of a photon increases the atomic kinetic energy by E r , so that the 
rate of increase of the energy is E = ^ E r = &b x 5 nK/s. Consider now the case of 40 K as 
an example of fermionic species. The fine structure is significantly reduced with respect 
to rubidium, leading to a lower merit factor, M = 2 x 10 5 . Choosing again ft e ff = c j r , 
the scaling given above leads to a scattering rate that is 10 times larger than for Rb 
atoms, and the heating rate reaches E = x 120 nK/s. This may be too large for a 
reliable production of strongly correlated states based on Raman coupling. Note that for 
a precise comparison between the heating rates for Rb and K, one should also specify 
the polarisations of the Raman beams, and take into account the difference between 
the Clebsch-Gordan coefficients originating from the difference in the nuclear angular 
momenta [25]. 


Lanthanides: erbium, dysprosium. Consider now atomic species like erbium or dyspro¬ 
sium, which have a more complex electronic structure than alkali-metal species. These 
atoms have two outer electrons and an incomplete inner shell (6s 2 and 4f 10 for Dy). 
Because of this inner shell, the electronic ground state has a non-zero orbital angular 
momentum (L = 6 for Dy). The lower part of the atomic spectrum contains lines corre¬ 
sponding either to the excitation of one of the outer electrons, or of one electron of the 
inner shell. 
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Fig. 6. - A few relevant levels of Dysprosium. The blue resonance line is split into three 
components by fine structure coupling. In order to minimize heating due to the random recoils 
associated with the spontaneous emission of photons, a Raman coupling between the Zeeman 
sublevels of the electronic ground state can be generated by a laser close to resonance with the 
narrow red transition at 741 nm. 


We show in fig. 6 a few relevant levels for dysprosium (for erbium, see [26]). The 
(broad) blue line corresponds to the resonance transition 6s 2 —6s 6p with the natural 
width T^/ 2,tt = 30 MHz; it is split in three components by fine structure coupling. There 
are also several narrow transitions such as the red line at 741 nm shown in fig. 6, cor¬ 
responding to the excitation of an electron of the inner shell (4f 10 —> df^d, term 5 K°, 
width T r / 27 T ~ 2 kHz) [24]. We irradiate the atom with a monochromatic laser in order 
to induce a Raman transition between ground state sublevels. This laser is characterized 
by its detunings A^ r from the blue and red transitions, and the corresponding Rabi 
frequencies /^ ;r . As we show below, one has to choose a small red detuning |A r | <C |A&| 
in order to maximize the merit factor. 

For the relevant choice of parameters, the dominant term in the Raman coupling 

2 

comes from the red transition, ~ > an d the spontaneous emission rate is obtained 

by summing the contributions of the blue and red lines( 3 ): 


(67) 



+ r 6 



The Rabi frequencies Kb jV associated with the red and blue transitions are proportional 


( 3 ) In order to simplify the algebraic expressions, we give here results obtained within the 
rotating wave approximation, which is only marginally correct for this large detuning with 
respect to the blue resonance line. 
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to the reduced dipoles d of these transitions, which are themselves proportional to 
y / T^ 7 , hence: 


( 68 ) 



V 

T b ' 


Using the fact that the ratio Tr/T^ is small (~ 10 4 ), an approximate calculation leads 
to the maximal merit factor 


(69) 


M ~ 


IV 


which is obtained when the detunings are such that 


(70) 


I Ar | 

I a 6 | 



The optimal detuning of the Raman coupling laser is only a fraction of angstrom from 
the red line, while A 5 is of the order of an optical frequency. This leads to a situation 
that is much more favorable than for alkali-metal species [compare eqs. ( 66 ) and (69)], 
with a merit factor A4 ~ 10 7 and a scattering rate 7 ~ 10 -3 s -1 for a Raman coupling 
equal to the recoil energy. The residual heating is then E ~ k b x O.lnK/s for both 
bosonic and fermionic isotopes. 


3. — Non-Abelian potentials and spin-orbit coupling 

We now generalize the problem of simulating external magnetic fields and ask if it 
is possible to take advantage of the internal structure of an atom to generate a single¬ 
particle Hamiltonian like 


(71) 


H = 


(p ~ M)) 

2 M 


where A(r) is not anymore just a vector field, but an ensemble of three matrices A = 
(A x , Ay, Az) acting in the internal Hilbert space of the atom: 

(72) A x (r) = [4 m ' n V)l, Mr) = K m ’ n) (r)], A z (r) = [M’ n) (r)}, 


where the indices m, n run over a basis of this Hilbert space. 

The non-Abelian character is linked to the fact that in a given point r, two components 
generally do not commute: 


( 73 ) 


[A x {r),A y (r)\ ^ 0 . 
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It is clear that this non-Abelian character can enrich the physics of gauge fields. Con¬ 
sider for example a Aharonov-Bohm experiment; in the Abelian case, the presence of 
a magnetic field entails that the state of the particle acquires a geometric phase when 
the particle travels around a closed contour C. In the non-Abelian case, this phase is 
replaced by a matrix acting in the internal Hilbert space: even if the atom moves arbi¬ 
trarily slowly, it may not end up in the same internal state after traveling around the 
contour ( 4 ). 

3*1. Non-Abelian potentials in quantum optics. - The idea to generate non-Abelian 
potentials in a geometric manner was proposed by Wilczek and Zee [27]. In this work the 
authors considered a generalization of the adiabatic theorem in the case where the Hamil¬ 
tonian possesses a group of eigenstates that remain degenerate (or quasi-degenerate) in 
the course of the evolution, this group being well-separated from the other eigenstates. 
This analysis triggered many studies ranging from molecular physics to condensed matter 
physics [28, 29]. Here we will focus on atomic physics implementations, initially proposed 
in [30, 31, 32, 33] (for reviews, see [4, 2]). 

In the quantum optics context, we suppose that the internal atomic Hilbert space is 
a subspace £ q of dimension q that is well separated from the rest of the spectrum. We 
denote {|?/q(r)), • • • ? |V^( r ))} an orthonormal basis of £ q , composed of eigenvectors of 
the internal atomic Hamiltonian. The adiabatic approximation assumes that the state 
of the atom can be written at any time as 

(74) ^(r) « y^0„(r) \ip n (r)). 

n— 1 

We project the full Schrodinger equation on the subspace £ q and obtain q coupled equa¬ 
tions: 


(75) 


at 


(P - Mr )) 2 

2 M 


+ E(f) + V(f) 


®(r,t), 




\M r R) 


Here E stands for the q x q diagonal matrix formed with the energies E n (r ) and V(f) 
is a q x q matrix representing a geometric scalar potential, generalizing eq. (46). The 
most important term for our discussion is Berry’s connection *A(f), which is also a qx q 
matrix with components: 

(76) = (^ n ( r )|vV> m (r)), 


( 4 ) The precise definition of a non-Abelian gauge field considers two contours C\ and C 2 and 
states that the final internal state of the particle should not be the same when the particle 
travels over C\ and then C 2 , or the reverse [2]. 
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Fig. 7. - (a) Tripod configuration leading to the emergence of a non-Abelian gauge potential. 
An atom has its electronic ground state with angular momentum J g — 1 coupled to an electronic 
excited state with angular momentum J e = 0. (b) Possible laser configuration with three beams 
(plane waves) at 120 degrees, with orthogonal linear polarizations. Each beam drives one of the 
three transitions gj H e. (c) Dressed state picture: one linear combination of the gj sublevels 
(bright state) is coupled to e. The orthogonal subspace of dimension 2 forms a “dark” subspace, 
in which the dynamics is governed by a non-Abelian gauge field. 


3*2. Tripod configuration and 2D spin-orbit coupling. - As a concrete example consider 
the tripod geometry represented in fig. 7, where a ground state with angular momentum 
J g = 1 is coupled to an excited state with angular momentum J e = 0. The coupling is 
provided by three laser beams with wave vectors kj (j = 1, 2, 3), with the same frequency 
and with orthogonal linear polarizations, propagating in the xy plane at 120 degrees from 
each other. By construction only one linear combination of the ground state manifold is 
coupled to the excited state. This bright state reads: 

(77) |B(r)> = T (g-ifei-r^i) + + e-^lfla)) , 

where | gf) is an orthonormal basis of the J g = 1 manifold, chosen such that the laser 
beam j drives the transition gj o e. The energy of the bright state is shifted from 
its initial position by the coupling with the light. The orthogonal dark subspace of the 
ground state manifold forms the desired space £ q , with here q = 2. After some algebra 
one reaches the expression of the non-Abelian Berry’s connection associated to this dark 
subspace [34]: 


(78) 


hk , „ „ v 

2 \Px'U'x T Oy'U'y) 



where the Oj (j = x, y , z) stand for the Pauli matrices and where u x , u y are orthonormal 
vectors in the xy plane containing the three wave vectors kj . We then obtain the following 
Hamiltonian for an atom moving in the xy plane when its internal state lies in this dark 
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subspace: 


(79) 


H = 


(P ~ A(r)) 
2 M 


P 2 

2M 


hk 
2 M 


(PxVx +Py°y) + • • • • 


The Hamiltonian (79) constitutes a prototype of spin-orbit coupling (SOC). The or¬ 
bital motion of the atom, here its 2D linear momentum p x , y , is coupled to its (pseudo) 
spin degree of freedom described by o x ,y We note that this SOC originates from the re¬ 
coil/Doppler effect; this is very different from the standard SOC in atomic or condensed 
matter physics, which has a relativistic origin. In standard SOC a charged particle (elec¬ 
tron) moves with velocity v in a region with an external electric field £. In the reference 
frame of the particle, a motional magnetic field & = v x £ appears; SOC then results from 
the interaction between this motional magnetic field and the intrinsic magnetic moment 
of the particle, proportional to its spin, p = 7 S'. In atomic physics, the relevant electric 
field is the Coulomb field generated by the nucleus and it has a radial structure £ || r. 
This leads to a SOC term oc (rxp)-S = L-S. In solid state physics, the electric field can 
be considered as uniform and leads to the so-called Rashba or Dresselhaus Hamiltonians, 
depending on the origin of £ (for a review, see [3]). The coupling generated in eq. (79) 
simulates such a 2D Rashba-Dresselhaus Hamiltonian. 

In solid state physics, SOC is at the origin of many applications, in particular in 
relation with spintronics [3]. It also plays an important role in fundamental physics 
since it is at the origin of the notion of topological insulators, which are insulating 
materials in their bulk but can conduct electricity at their surface [35, 36]. This surface 
electricity conduction is topologically robust, as for the Quantum Hall effect. However 
in contrast with the Quantum Hall effect, there is no external magnetic field that would 
break time-reversal invariance. In a Fermi gas with point-like interactions, it can lead 
to the formation of novel superfluids [37, 38] and possibly the creation of Majorana 
quasiparticles [39, 40, 41]. Finally we note that even at the single-particle level, the SOC 
Hamiltonian given in eq. (79) can have spectacular consequences, since it changes the 
nature of the spectrum. In particular the ground state, which is non-degenerate (p = 0) 
in the absence of atom-laser interaction, becomes massively degenerate since all states 
with momentum p with \p\ = hk/ 2 can be associated to a spin state such that they have 
a zero-energy. This massive degeneracy of the single-particle ground state changes the 
nature of Bose-Einstein condensation for an ideal gas [42]. 

3*3. ID version of spin-orbit coupling. - The multi-dimensional implementation of 
spin-orbit coupling has not yet (to our knowledge) been implemented experimentally. 
However a ID version of it, initially identified in [43], has been realized by several groups, 
first with bosons [44] and then with fermions [45, 46] . 

To describe this implementation, we consider the A scheme of fig. 3b and we suppose 
that the transitions gj e are driven by two plane waves with wave vectors ±k along 
a given direction, say x. We suppose that the laser excitation is far detuned from the 
ground to excited state resonance, and that we can eliminate the excited state in a 



24 


J. Dalibard 



Fig. 8. - Energy levels of the ID SOC Hamiltonian given in eq. (81), labeled by the momentum 
index p of the family (80). The detuning A is chosen equal to zero. The Rabi frequency is given 
by hKo/E r = 0 (dash-dotted), 1 (continuous), 4 (dashed) and 6 (dotted). 


perturbative manner. We then deal only with the internal subspace associated to the 
ground states g\ and $ 2 - We introduce the family of states labeled by the momentum p: 


(80) 


F(p) = {\gi,p-hk), \g 2 ,p + hk)}, 


which is globally invariant by the action of the atom-laser coupling. This coupling is 
described in this family by a 2 x 2 matrix: 


(81) 


H(p) 


(p-hk) 2 /2M- 

Hkq/2 


hA / 2 Hk o / 2 

(p + hk) 2 /2M — hA/2 


which can also be written in terms of the Pauli matrices: 

/ x tt/ \ 1 / i\ 2 hA ^ hK 0 ^ , j ^ 

(82) H(p) = — [p-Aj + -Y a z + — with A=hka z . 

This corresponds to the desired form, since the “vector potential” is an operator with 
respect to the internal degrees of freedom. 

It is clear that the three components of A commute with each other in this particular 
case, so that it does not correspond to a non-Abelian coupling. However in spite of 
this restriction, some characteristics of SOC remain relevant. For example, for A = 0 
and a relatively small amplitude of the Raman coupling, one still obtains a non-unique 
single-particle ground sate. We plotted in Figure 8 the energy levels obtained from the 
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diagonalization of the Hamiltonian (81): 

(83) E ± ( P ) = £ i+ E r ±l 

The lowest level, E-(p), has two symmetric minima in p and — p if the Rabi frequency 
Ko is such that 

(84) Hkq < 4 E r . 

Above this value, the minimum of E_ (p) is located in p = 0. In the limit kq E T , 
we recover the situation where the adiabatic approximation is valid [cf. eq. (61)]. The 
interesting case for SOC, z.e., a degenerate ground state, thus corresponds to the regime 
opposite to the one studied in subsect. 2‘2. 

At the single-particle level, the achievement and the description of ID SOC is thus 
quite simple. It is essentially a reinterpretation of the usual dressed atom diagram, 
in which the coupling between internal and external degrees of freedom is due to the 
Doppler shift in a plane running wave. At the many-body level, SOC is at the origin 
of several interesting phenomena even this conceptually simple ID configuration (for a 
review, see [47]). Several questions are still open concerning the phases that may appear 
for bosons (see for example [48]) as well as for fermions, with the possibility to generate 
Majorana particles at the edges of a chain of atoms in a topological superconductor 
[49, 50, 51, 52, 53, 54, 55, 56, 57]. 



4. — Gauge fields on a lattice 

In the two previous sections we considered the generation of a gauge field for an atom 
moving freely in space. Another important class of problems deals with the magnetic 
phenomena that can appear in the presence of a spatially periodic potential. This ques¬ 
tion emerges in particular when one studies the effect of a magnetic field on the electron 
fluid of a crystal. 

The richness of the problem is linked to the existence of two length scales, which can 
“compete” which each other and lead to frustration phenomena. The first length scale 
is the period a of the potential (in the following we consider a 2D square lattice). The 
second length scale is the magnetic length that we already introduced ^ mag = \/h / qB. 
The ratio of these two length scales can be written in terms of the flux <f> = Ba 2 of the 
magnetic field through the unit cell of the lattice and of the flux quantum 4>o = h/q\ 


(85) 


a 2 _ qBa 2 _ o 4> 


The ratio a 2 /^ ag is thus equal to the Aharonov-Bohm phase accumulated by the particle 
when it travels along the sides of the unit cell. 
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As long as the flux <h remains small compared to <f> 0 , one does not expect spectacular 
modifications compared to the free particle case: the lattice step is much smaller than the 
typical size of a cyclotron orbit, which is thus only marginally affected by the discretiza¬ 
tion of space. This regime 4> <C 4>o corresponds to the case of “ordinary solids” (lattice 
period of a few Angstroms) and of magnetic fields that can be achieved in the laboratory 
(a few tens of Teslas). On the opposite, for large magnetic fields such that 4> ~ <hch the 
competition between these two length scales is at the origin of new phenomena, such as 
a fractal structure for the single-particle spectrum. This regime can be explored with 
synthetic materials (see e.g. [58, 59, 60]) or with artificial magnetic fields, in particular 
with cold atoms in an optical lattice. 

A key notion that will appear in the following is the topology of an energy band. 
It is characterized by an integer number, the Chern index, which indicates if a given 
filled band can contribute to a particle current along a given direction, say y , when a 
force is applied along the orthogonal direction x. This notion is directly inspired from 
the study of Hall conductivity [61]; we will introduce it here by considering the case of a 
square lattice placed in a uniform magnetic field, the so-called Hofstadter model [62], that 
was recently implemented with cold atoms in optical lattices [63, 64]. This topological 
characterization of bands via their Chern index can be generalized to other classes of 
lattice configurations, such as the honeycomb model introduced by Haldane [65] that has 
also been realized recently with cold atomic gases [66]. 

4T. Tight-binding model. - In most of this section, we will consider for simplicity a 
single-band, tight-binding model in a two-dimensional square lattice (Hubbard model). 
We denote by |j, l) the state of the particle localized on the site r = a(ju x +lu y ) (j, l G Z) 
and we assume that the particle moves in the lattice via tunneling from a given site (j, l) 
to the four neighboring sites j =b 1, l and j, l ± 1. Denoting the tunnel amplitude by J, 
the Hamiltonian in the absence of magnetic field is: 

(86) h = -J^2(\j + M)0AI + LM + i)0AI) + h.c. 

3d 


The eigenstates of this Hamiltonian are the Bloch states \ip(q)) labelled by the Bloch 
vector q = (q x ,q y ) with energy E(q) 

(87) | ip(q)) = 'Y^e ia ( 3 qx+lq v )) \j,l), E(q) = -2 J [cos (aq x ) + cos(a^)]. 


The energy spectrum is thus a band centered on E = 0 with a full width 8J. Since 
two Bloch vectors differing by a vector of the reciprocal lattice correspond to the same 
eigenstate, we can restrict the domain accessible to q to the first Brillouin zone —it/ a < 
qj < 7r/a, j = £,;?/, in which case all | 'ip(q)} are independent and form a basis of the 
Hilbert space. 
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We take into account the presence of the magnetic field by assigning a complex value 
to the tunneling matrix elements: 

(88) -J\j + l,l)(j,l\ —► >J~U) \j + l,l},(j,l\ 

and similarly for the links in the other direction. We also set (j){j + 1,/ j,l) = 
j0 1,/) to ensure that the Hamiltonian remains Hermitian. This technique 
of using complex tunnel coefficients is well known in condensed-matter physics and is 
called the Peierls substitution. In the case of a real magnetic field B(r ), the Peierls 
substitution states that the phase 0(j, l —> j + 1, l) is calculated (once a gauge choice has 
been made) using the vector potential: 


(89) (j)(r r') = ^ J A dr. 

The justification of the Peierls justification is far from trivial and we refer the reader to 
[67, 68] for a discussion. However in the context of cold atoms, the possible difficulties 
associated to this justification are not relevant; indeed one looks for a direct implemen¬ 
tation of the complex tunnel matrix elements of eq. (88), without deriving them from a 
continuous model with external B(r) and A(r). 

As for the continuous case, there is a gauge freedom for the spatially discretized 
problem. An infinite number of choices for the phases {4>(j, l —> j +1, /), (j){j, / —> j, l +1)} 
lead to the same physical situation, and one can switch from one set to the other by 
changing the phases of the | j, /)’s. One can show quite generally that the gauge invariant 
quantities for the discrete problem are the sums: 

= j + 1,1) + 4>{j + 1, l j + 1, l + 1) 

(90) + (/>(j + ld + l^ j, Id- 1) + l + 1 j, l). 

Each sum represents the accumulated phase when the particle travels counterclockwise 
around the cell having the site (j, /) in its lower left corner. Using eq. (89), one sees that 
is just the flux of the magnetic field B(r) through the considered plaquette. 

We will be interested in the following in the case of a uniform field on the lattice, be., 
the same magnetic flux <f> across each lattice cell. Using eq. (89) with the Landau gauge 
A = — Byu x , we find that this amounts to take 

(91) l j,l + 1) = 0, l j + 1, l) = —2na l. 

One can immediately check that eq. (90) yields in this case the same flux 
for all cells. With this gauge choice, the target Hamiltonian reads 

R=-J^(e- i2 ™ i |j + 1,0001 + \j,l + l)(j,l\) + h.c. 

3,1 


( 92 ) 
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4*2. Hofstadter butterfly. - Before looking at the possible implementations of the 
Hamiltonian (92) with a cold atom setup, we briefly discuss the main properties of its 
spectrum and of its eigenstates. First we note that this Hamiltonian is invariant in the 
change a —> a + 1. Therefore it is enough( 5 ) for this spatially discretized problem to 
study the spectrum of the Hamiltonian for a between 0 and 1 . 

For the gauge choice leading to the Hamiltonian (92), the presence of the magnetic 
field breaks the translational invariance along y , but preserves the invariance along x. 
Hence one can still look for the eigenstates of H as Bloch functions along x : 

(93) W = XT* e iia9 *|j,0> 

3d 


where the coefficients Ci are unknown at this stage. Inserting |T) in the eigenvalue 
equation for H, we find the recursion equation called the Harper equation 

(94) Ci— i + 2 Ci cos(27r<a l + aq x ) + C /+1 = ~(E/ J ) Ci. 

One is then left with the diagonalization of a tri-diagonal matrix. 

An important simplification occurs when a is a rational number, a = p'/p, where 
p,p' are coprime positive integers. Then one recovers a periodic problem also along p, 
but with an increased spatial period pa instead of a. Indeed one gets in this case 

4>(j, l + p -m j + 1, l + p) = —2na(l + p) = —27 ral — 2irp' 

(95) = (f)(j , l —► j + 1, Z) modulo 27 t. 

We are then back to the problem of the motion of a particle on a periodic lattice, but 
with a unit cell of size a x (pa) and p non-equivalent sites per unit cell. In this case the 
general result is that the initial energy band (of width 8 J) gets fragmented in p subbands, 
generally separated by forbidden gaps( 6 ). 

We choose the case a = 1/3 asa concrete example and refer the reader to Appendix 2 
for a general discussion. The unit cell has a length a along x and 3a along y. It contains 
three sites denoted \A),\B),\C) in fig. 9. Each cell is labelled by the indices j, where 
we take by convention the A site at point a [ju x + 3l'u y \. We look for the eigenstates of 
the Hamiltonian as Bloch functions labelled by the Bloch vector q = (q xi q y )\ 

(96) |*(g)> = e i a(jq x +iVq v ) (^1 A u ,) + ,,) + e 2ia ^/3 3 \C jt ,,» 

3,1' 


( 5 ) One can even restrict the domain of study to the interval a G [0,1/2], since one switches 
between a and -a by changing the orientation of the x axis. 

( 6 ) These gaps may be replaced by a contact with a Dirac point, like around Em 0 for a even. 
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Fig. 9. - Left: Choice of a unit cell of size a x 3a for a magnetic flux a = 1/3. Right: the three 
subbands resulting from the fragmentation of the band of width 8 J. 


with \A jX ) = \j,3l'), | B j<v ) = \j,3l' + 1), | C jX ) = \j,3l' + 2) and 

(97) q x e] - Tr/a,n/a], q y e\ - n/(3a),n/(3a)]. 

The coefficients /3 r are obtained from the eigenvalue equation H |\L) = E |T), which 
amounts to look for the eigenstates of the Hamiltonian in reciprocal space: 

/2cos (aq x ) e [aq y e~ [aq y \ 

(98) 77(g) = — J | e _ia<3 ^ 2 cos(aq x + 2tt/ 3) e iaqy I 

V e iaqy e~ iaq y 2 cos(aq x + 4tt/3)/ 

written here in the basis {\A),\B),\C)}. The diagonalization of this 3x3 matrix gives 
three eigenvalues, which are functions of q. When q x and q y vary in the intervals (97), 
one finds the three energy bands indicated in fig. 9. The energy spectrum is symmetric 
with respect to E = 0 and its total width is reduced with respect to the value 8 J found 
for a = 0. The lowest subband has a width ^ 0.7 J and it is separated from the middle 
band by a gap of ~ 1.3 J. 

The spectrum of H for arbitrary a is represented in Figure 10. It has a very specific 
fractal structure called Hofstadter butterfly [62]. One can understand the origin of this 
structure by comparing two close values of the flux like 1/3 and 10/31, and noting that 
they lead to very different results, the first one with three subbands (fig. 9), the second 
one with 31 subbands. 

It is interesting to connect the result obtained in the lattice case for low magnetic 
fluxes (a 1) to the Landau level structure of the continuum case. A zoom on the lower 
left corner of the Hofstadter butterfly is shown in fig. 11. One sees that each bandwidth 
gets very small compared to the band spacing, so that we recover in the limit a 0 
the notion of well defined energy levels. The levels are approximately equidistant with 
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a 


Fig. 10. - Hofstadter butterfly: energy spectrum for a particle moving on a square lattice in 
the tight-binding approximation in the presence of a uniform magnetic field. The flux <f> of the 
magnetic field across a unit cell is such that <f> = a<f>o, where <f>o = h/q is the flux quantum 
[62]. The calculation has been done for rational values of a, a = p /p, with denominators p up 
to 100. 


an energy that can be written E n « —4J + (n + 1/2) Hlj c , with uj c = qB/M e ff. Here the 
effective mass M e q is obtained from the dispersion relation at the bottom of the band 
[see eq. (87)]: M e ff = h 2 /(2Ja 2 ). 

4*3. Chern number for an energy band. - The subbands that result from the frag¬ 
mentation of the initial band in the presence of a magnetic field have specific, non-trivial 
topological properties. These properties are at the origin of a quantized Hall conductance; 
if one applies a constant force F along one axis of the lattice, say x, then a current of 
particles flows in the orthogonal direction. Suppose for example that the lowest subband 
is filled with fermionic particles, while the others subbands are empty. In the absence of 
topological properties, this would correspond to a band insulator and no current would 
be expected. However due to the topology of the lattice, a quantized non-zero current 
along y appears, which is characterized by the so-called Chern index C, an integer number 
associated to the considered subband. The Chern index is defined such that C particles 
cross a given x-oriented link of the lattice during the time interval tp = h/aF. This time 
£b is the period of the Bloch oscillation phenomenon that occurs in the presence of the 
force F for the lattice of spatial period a. The proof of this remarkably simple result, 
as well as relevant references, is given in Appendix 2. The important relation between 
the topological bulk property described by the Chern index and the edge currents that 
appear in a finite-size sample is addressed in [69, 70]. 
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Fig. 11. - The Landau levels recovered for the Hofstadter butterfly (zoom on the lower left corner 
of fig. 10). The low energy levels form a quasi-discrete spectrum with E n /J —4+47r(n+l/2)o;, 
n integer, corresponding to the Landau spectrum in (n + l/2)huj c , where the cyclotron frequency 
uj c is determined using the effective mass M e ff = ti 2 /(2Ja 2 ). 


A Chern index can also be associated to the Landau levels of a free particle in the 
presence of a uniform field (subsect. 1‘2). In particular the Chern index of the lowest 
Landau level is C = 1. Consequently the lowest subband in a lattice configuration is said 
to be topologically equivalent to the LLL if it also has C — 1. This is the case of the 
lowest subband of the Hofstadter butterfly when one chooses a = 1/p with p integer. 

As explained in Appendix 2, the Chern number is equal to the integral over the 
magnetic Brillouin zone of the Berry curvature B(q) associated to the variations of the 
Bloch states \^f q ) with the Bloch vector q. The Berry curvature can itself be determined 
from the dynamics of an atomic wave packet in the lattice, as proposed in [71, 72, 73]. 
This technique was recently implemented by the Munich group (see the contribution of 
Immanuel Bloch to this volume and [74]) using a topological band produced by lattice 
modulation, which we describe in the next section. 

5. — Generation of lattice gauge fields via shaking or modulation 

We now turn to the description of some possible ways to generate non-real tunneling 
matrix elements for an atom moving in an optical lattice. Many possible implementations 
have been proposed in the literature and our presentation will not be exhaustive. In 
particular we will not describe methods based on lattice rotation, for which the basic 
principles are similar to the bulk case described in subsect. 1*4 (see e.g. [75, 76, 77, 78] 
and refs. in). 
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5*1. Rapid shaking of a lattice. - We consider in this section a Hamiltonian which 
depends explicitly on time and corresponds to a shaken lattice. The potential reads in 
the one-dimensional case 

(99) V(x,t) = V[x - x 0 (t)], 

where xo(t) is a time-periodic function with period T. We will suppose in this subsection 
that T is very small compared to the other time scales of the problem, such as h/J. A 
rigorous treatment of the problem is based on a Floquet-type approach. More precisely, 
one can develop a method that separates in a non-ambiguous way the dynamics related 
to an effective, time-independent Hamiltonian and the micro-motion caused by the fast 
modulation xo(t) [79, 80, 81]. Here for simplicity we will not follow such a rigorous 
treatment, but we will take a time-average of the Hamiltonian H(t), once it has been 
written in a suitable form. The link between this simple-minded approach and rigorous 
treatments is made in [80], for example. 

It is not immediate to transcribe the Hamiltonian of a particle in a shaken lattice 

^2 

(100) H 1 (t)*tJL- + v[x-x 0 (t)] 

in the tight-binding approach. The simplest strategy is to use the unitary transformation 
generated by 

(101) U(t) = exp[ix 0 (t)p/h], 
so that the Hamiltonian after transformation 


A(t ) = Mx 0 (t ). 

This amounts to applying a time-dependent vector potential in the direction of the shak¬ 
ing^). Then the tight-binding version of the Hamiltonian is obtained in a straightforward 
way: 

(104) H(t) = — J e [Ma± o(t)/h 5> + l>0l+h.c. 

3 


( 102 ) 

reads 

(103) 


H{i) = U(t) fh(t)U ] (t) + 


^ = —2 


( 7 ) The term A 2 (t)/2M is a mere time-dependent global shift of the energies. It can be elimi¬ 
nated by a unitary transform and has no physical consequence. 
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Now using the Magnus expansion at the lowest order (see e.g. [80] and refs, in), we 
simply time-average this Hamiltonian over one period of the fast oscillation. We obtain 
the renormalized tunnel matrix element: 


(105) 


J — J (^Maxo{t)/h\j 


In the context of cold atom physics, this result has first been considered as a way to 
change the magnitude of the tunnel matrix element by using a sinusoidal variation of the 
lattice displacement xo(t) [82]: 

(106) ^-x 0 (t) = £ 0 sin (fit + 0) J = J Jo(£ 0 ), 

where J n designates the n-th Bessel function of the first kind. This modification of J in 
a ID vibrating lattice has been accurately checked by the Pisa group [83]. 

Here we are interested in obtaining complex tunnel matrix elements, which is not the 
case for the sinusoidal modulation of eq. (106). A non-symmetric temporal modulation 
of xo(t) provides the desired result [84]. For example one can divide the time-period T 
in two unequal time T\ and X 2 (Xi + X 2 = T) and have the lattice moving at uniform 
velocity Vi during Ti (i = 1,2) with v\T\ + V 2 T 2 = 0 so the displacement of the lattice 
is zero in one full period T. Then the renormalized tunnel coefficient obtained from this 
saw-tooth modulation reads 


(107) 


J ^\ Maxo(t)/h ^ ^iMavi /h _|_ Mav 2 /h 

J T T 


which is generally non-real. This technique was implemented by the Hamburg group in a 
ID lattice, and the displacement in momentum space of the bottom of the lowest band, 
which constitutes a signature of a non-real effective tunnel coupling, was observed [84]. 
A different implementation of a ID tight-binding Hamiltonian with complex tunneling 
coefficients, based on a combination of Raman coupling and radio-frequency fields, was 
realized at NIST [85]. 

The technique of a fast global shaking of the lattice can be generalized in two dimen¬ 
sions to generate a non-zero flux through the cells of the lattice. A close look at the 
effect of the modulation shows that it does not lead to the desired flux if the sides of the 
unit cell are parallel to each other. In this case, the phases on two parallel sides indeed 
cancel each other, when one looks at the total phase accumulated over the border of a 
unit cell as in (90). On the contrary, this method does provide a non-zero flux for a 
triangular lattice, as demonstrated by the Hamburg group [86]. Note however that the 
obtained flux is not uniform but staggered, with equal magnitude and opposite sign for 
two adjacent triangles. Using also a fast modulation and starting from a honey-comb 
type lattice, the ETH group could recently implement experimentally the Haldane model 
[65, 66]. 
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5‘2. Resonant shaking/modulation. - Another possibility for generating complex tun¬ 
neling matrix elements is to shake the lattice system at a resonant frequency. To reach 
this goal, one first adds to the lattice an external, static potential that lifts the degeneracy 
between adjacent sites by a energy offset Ml. Then one shakes the lattice at frequency £2, 
and the phase of the function describing the shaking gets “printed” on the corresponding 
tunnel matrix element. To prove this result we follow again a simple approach based 
on the lowest order of the Magnus expansion. A more rigorous, systematic expansion in 
powers of 1 /Q can be found in [87, 81]. 

Let us consider for example the case of a ID lattice, where the energy offset Ml between 
adjacent sites is obtain by applying a linear potential — Mix/a , where x is the position 
operator and a the lattice constant. The global shaking of the lattice is described by the 
function xo(t) as in eq. (100) and we proceed with two successive unitary transformations. 
The first one uses the same operator as in eq. (101) and leads (up to additive constants) 
to 

(108) Hi = -j]T I j + 1)01 + h.c.) - hnYi j liXil- 

3 3 

The second unitary transform eliminates the linear potential and replaces it by an extra 
time-dependence of the tunnel coefficients [see eqs. (148) and (150) in Appendix 2 for 
details]: 

(109) H(t) = -J £li + 1)01 + >“• 

3 

We assume the same sinusoidal modulation of the lattice as in eq. (106) and take the 
time average of the Hamiltonian (109) over one period of modulation 27r/D. We obtain 
the renormalized tunneling coefficient 

(110) J=J (e iKo sin(n t+4>)-SU]} = JJ^) gi 

This result differs from the one obtained for rapid shaking [eq. (106)] in two ways. First 
the zeroth order Bessel function Jo is replaced the first order function J\ ; in particular no 
tunneling occurs in the absence of modulation (J7i(0) = 0), because of the energy offset 
between adjacent sites. Second and most importantly, the phase </> of the modulation is 
now “printed” on the renormalized tunnel coefficient, which thus becomes complex. 

For a ID lattice, this printed phase can be viewed as a mere gauge transform in which 
one redefines the relative phases of adjacent Wannier functions. The procedure has wider 
applications in 2D since it can lead to a non-zero, uniform phase on each lattice cell. A 
first generalization of this method to the two-dimensional case was proposed in [88], but 
it presented some drawbacks discussed in [89]. An alternative method that provides the 
desired uniform flux over the lattice was proposed [90] and implemented by the Munich 
and MIT groups [63, 64]. They relaxed the condition of a uniform shaking of the lattice 
sites and implemented a modulation that varies from site to site. This implementation 
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Fig. 12. - Left: Infinite ladder with two sides separated by the distance a. The ladder rungs 
are equidistant, separated by a distance b. The particle can jump from one site to the next by 
tunneling. A uniform magnetic held creates a flux <f> on each cell. Right: Simulation of this 
ladder with a double optical lattice. A laser beam propagating along the lattice direction induces 
transitions | g) CA |e), similar to the tunneling along the ladder rungs in the left picture. The 
laser phase 0 = ky varies linearly with the rung index j and gets “printed” on the corresponding 
matrix element. 


was made possible using the dipole potential created by an extra pair of laser beams, 
which created the necessary modulation of the lattice sites( 8 ). 

6. — Generation of lattice gauge fields via internal atomic transitions 

The possibility to use internal states in an optical lattice brings some new flexibility 
for the generation of complex tunnel coefficients in a lattice. The process in play is laser 
assisted tunneling , which we will first present on a simple, one-dimensional system, before 
switching to the case of an infinite two-dimensional lattice. 

6*1. Laser assisted tunneling in a ID ladder. - We are interested here in a ladder with 
two parallel sides as in fig. 12. This ladder is submitted to a uniform magnetic field and 
we suppose that it is well described by the same tight-binding approximation as above. 
The magnetic field is characterized by the flux across a unit cell. 

To simulate this ladder we use an atom with two internal states | g) and \e). A 
trapping potential created by an optical lattice along y localizes the atoms along the 


( 8 ) Since the modulation is induced by an auxiliary laser field, this scheme is sometimes called 
laser-induced tunneling. Although this terminology is certainly correct in the present case, 
we prefer to reserve it for schemes where a laser beam is used to simultaneously (i) induce a 
transition between internal atomic states, and (ii) induce a jump of the atomic center-of-mass 
between adjacent lattice sites, as described in subsect. 6'1. 
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rungs of the ladder. We assume that this potential is different for | g) and |e), so that 
each internal state is trapped on a given side of the ladder (| g) on the left, \e) on the 
right). We suppose that these two internal states are stable (no spontaneous emission) 
and we consider a laser excitation that allows one to induce a resonant transition between 
them. As explained in subsect. 2*3, this atomic level scheme can be obtained with species 
with two outer electrons, as well as with alkali metals, if one uses a Raman transition 
between sublevels of the atomic ground state. 

Using a single-band model, we denote \j) g and |j) e , j G Z, the spatial states associated 
to the two internal states | g) and \e). To simplify the algebra we suppose that the 
corresponding wavefunctions are identical, except for the translation over a: 

(111) Wj(r) = wq(x, y — jb), Wj(r) = w 0 (x - a,y - jb). 

In the absence of laser excitation, we then have two independent ID lattices with the 
Hamiltonian: 

(112) H = H g + H e , H a = -J y ^2\j + l) a a (j\ +h.c., a = g,e. 

iez 


In order to couple these lattices, we use a laser beam that is resonant with the transition 
| g) GG |e), with a wave vector k parallel to the y axis. The atom-light coupling reads: 

~ Jl K 

(113) l/= T e i ^|e)( 5 |+h.c. 

Let us rewrite this coupling in terms of the spatial states | j) a . Using the closure relation 
for internal and external states in the single band approximation: 


(114) i = |0>(5l® ^2\j)gg(j I + |e)(e|® ^ \j)e e(j\ 






we obtain 

(115) 


r XG/.A A/ • 


with 

(116) Vjj' = J w j ( r ) elky w j'( r ) d 2 r • 

In the tight-binding approximation the wave functions re|(r) and ( r ) are well localized 
and on can neglect their overlap as soon as j ^ j' . This allows us to simplify the 
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expression of Vjjf 


(117) 


Vjj' ~ Sj 


,j' y J W 0 (x-a,y- jb) e lky w 0 (x, y - jb) d 2 r 


A. / Jjkb 

U j,j' d X C 5 


where the tunnel coefficient J x is proportional to the Rabi frequency and to the overlap 
of the Wannier functions for the internal states | g) and |e): 

JlK C 

(118) J x = — J w 0 (x - a, y) e lky w 0 (x, y) d 2 r. 

For a function wo(x,y) that is even in the change y —>> —y, this coefficient J x is real: 

r 

(119) Jx = Y w °^ x - a,y ^ cos ^ ky ^ w o( x ^y ) d2r - 

The expression (117) has the desired structure. The tunnel matrix element is complex, 
with a phase that increases linearly with the index j. The total phase accumulated along 
the contour of a cell 


( 120 ) 


7 \ J y v 


t „i(i + l )kb 


is non-zero and it takes the same value (j + l)kb — jkb = kb for all cells. One simulates 
in this way a uniform magnetic field over the whole ladder. 

6*2. Lattice with artificial dimension. - One can notice that in the preceding descrip¬ 
tion of a two-leg ladder, nothing requires the length a to be non-zero. The ladder can 
have a dual nature: real along the y direction, and fictitious along x. The requirement for 
a non-zero a will come only later, when we try to increase the number of sites along the 
x direction. Starting from this remark, it was proposed in [91] to extend the preceding 
reasoning to the case of an atom with N internal states, | g n ), n = 1,..., TV, by choosing 
an atom-light coupling that induces the transitions 

T 

( 121 ) \9n,j) XJL + \9n+l,j). 


This allows one to simulate a ladder with N legs, each leg being associated to one internal 
state. We illustrate this in fig. 13 for the case TV = 3, taking for example the three Zeeman 
states m = 0, ±1 of an atomic ground state with angular momentum 1. 

6*3. Laser-induced tunneling in a 2D lattice. - We now extend the scheme outlined 
above to the two dimensional case. Now the length a is not fictitious and we consider 
a series of parallel ID lattices, corresponding alternatively to the internal states | g) and 
|e). One forms in this case a 2D lattice in which a site (j, l) is occupied by an atom in 
state | g) (resp. |e)) if j is even (odd). 
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Fig. 13. - Simulation of a double ladder for an atom with three internal states [91]. This scheme 
can give rise to a single-particle spectrum similar to the Hofstadter butterfly and to edge states 
with opposite currents along the “edges” associated to the states \gi) and ^ 3 ). 


The analysis of the phases associated to laser induced tunneling shows that the lattice 
then corresponds to a staggered flux. We recall that a uniform flux is obtained by taking 

(122) |j,Z)°-^li + M> 

whereas we have here 

i Ikb 

(123) | j, l) —» \j + 1, l ) if j even 
corresponding to the transition | g) —> |e) and 

-ilkb 

(124) | j,l) —» \j + 1, l) if j odd 

corresponding to the transition \e) —> \g). The sign of the accumulated phase on the 
contour of a cell alternates from one column to the next. 

This lattice with a staggered flux has interesting properties when interactions are 
taken into account [92], but it does not correspond to the desired simulation of uniform 
magnetism with non-trivial topological properties. One thus needs to find a way to 
“rectify” the magnetic field, in order to obtain a flux with the same sign over each cell. 
Several techniques have been proposed to reach this goal. The first one [93] consists in 
adding an extra linear potential in order to lift the degeneracy between the transitions 
| g) —> \e) going to the right (2 j —)> 2 j + 1) and those going to the left (2 j —» 2 j — 1). 
Another version, well adapted to alkaline-earth species, uses a super-lattice along the x 
direction with a spatial period 2 a [94]. In both versions, the resulting different transitions 
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Fig. 14. - Realization of a square lattice with a staggered flux, using a single light beam to 
induce laser-assisted tunneling. 


| g) —> \e) are provided by laser beams propagating along +u y and — u y , ensuring that 
the phases accumulated on all cells of the lattice always have the same sign. 

6*4. Optical flux lattices. - We now present a different way to simulate magnetism 
on a lattice, which consists in imposing that the light configuration is fully periodic and 
stationary in time. We consider a situation where all laser beams are monochromatic 
and form standing waves, and we ask if some interesting properties, i.e. energy bands 
with non-trivial topological features, can be found in this case. We will see that this is 
indeed the case, as it was first proposed by N.R. Cooper in 2011 [95]. 

We note first that due to the periodicity of the lattice, the phase accumulated by a 
particle when it travels along the sides of a unit cell must be a multiple of 27r. In the 
tight-binding limit, we know from our study of subsect. 4’2 that for such a phase, the 
physics is the same as for a zero magnetic flux, hence it is not relevant for our present 
goal. From now on, we will therefore assume that the potential created by the laser 
light is relatively weak, so that the lattice operates outside the tight-binding regime. 
In practice this means that the Rabi frequency characterizing the atom laser coupling 
should be on the order of (or smaller than) the recoil frequency E r /h. 

To proceed, we consider a situation where the internal dynamics is well described by 
the two-level approach of subsect. 2‘3, and we suppose that the parameters f2, 0 and e 1 ^ 
entering in the 2x2 matrix for the atom-laser coupling are periodic functions of space. 
Assuming for the moment that the adiabatic approximation holds, the vector potential 
A(r) -Berry’s connection- given in eq. (52) is also a periodic function of space. The 
contour integral of Aflr) over the unit cell must then be zero. If the Stokes theorem 
holds, the flux of Berry’s curvature B(r) through the unit cell is also zero, meaning 
that the lattice does not create any average magnetism. However there exist situations 
where the Stokes theorem is not valid, leading to a non-zero magnetic flux in spite of the 
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periodicity of A: this is precisely the definition of an optical flux lattice( 9 ). 

To obtain a “violation” of Stokes theorem, one needs to have singularities of A(r), 
ensuring that 


(125) 



B z dxdy. 


Let us take an example extracted from [95], in which the coefficients of the 2x2 matrix 
giving the atom-laser coupling vary as 

(126) Cl cos 0 oc sin(kx) sin(ky), D sin# e 1 ^ oc cos(kx) + icos(%), Q > 0. 

The unit cell of the lattice is —it < kx , ky < i r, and the phase 0 is ill-defined at the 
four points in the cell where sin# vanishes: ( kx,ky ) = (±7t/2, ±7t/2). This entails that 
Berry’s connection A, whose expression contains the term ±(h/2) cos# V0 [see eqs. (52) 
and (56)], may also be singular in these points. 

There are various ways to handle these singularities and determine the magnetic 
properties of this lattice configuration. One option is to evaluate directly B from eq. (53), 
and check that B is actually non singular at those points. A second possibility is to 
make different gauge choices A 1 , A 11 ,... over different regions of the unit cell, such that 
no singularity occurs over any region. The gauge transformations at the boundaries 
between these regions provide the desired flux( 10 ) [96]. A third approach consists in 
taking a contour for A that avoids the (discrete) singularities of </>, so that the equality 
in eq. (125) holds. The flux of B is then given by the contribution of these singularities, 
each being a multiple of irh. Still another method is to use the concept of Dirac strings 
to account for these singularities [97]. All these approaches of course lead to the same 
result, with a flux of B that can be non-zero and that is always a multiple of 2irh. 

The quantitative optimization of an optical flux lattice configuration amounts to make 
it as similar as possible to the case of a particle moving freely in a uniform magnetic 
field, in which case the energy eigenstates group together to form massively degenerate 
Landau levels (sect. 2). For this optimization, it is preferable not to rely on the adia¬ 
batic approximation since the relatively small laser coupling makes the validity of this 
approximation marginal [see eq. (61)]. Consequently, one solves the band problem for the 
Hamiltonian (40) and tries to fulfill simultaneously three conditions, focusing for example 
on the lowest band: (i) The band should have a non-trivial topology, characterized by a 


( 9 ) There is a duality between the concept of an optical flux lattice in real space, and the concept 
of a Chern band in momentum space. In the case of a Chern band, the physical problem is 
periodic over the Brillouin zone but one can still have a non-zero flux of Berry’s curvature B(q), 
corresponding to a non-zero Chern index C [cf. eqs. (167) and (171) in Appendix 2]. 

( 10 ) For the coupling (126), the gauge choice (52) leads to no singularity for A at the two 
points (kx, ky) = (7t/2,7t/2) and (— ir/2, — 7t/2) because (cos# — 1) vanishes in these two points. 
Similarly the gauge choice (56) leads to no singularity for A in the two other points (tt/2, —7t/2) 
and (— 7t/2,7t/2) because (cos# + 1) vanishes in those points. 
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e/e r 

Fig. 15. - Optical flux lattice configuration for an atom with a spin 1/2 ground state. Top left: 
A magnetic field lifts the degeneracy between the two Zeeman sublevels g±. Top right: Three 
identical laser beams propagate in the xy plane at 120 degrees from each other, with a linear 
polarization at an angle 0 with respect to the z direction. A fourth laser beam propagates along 
the z axis and induces, together with the three horizontal beams, a resonant Raman coupling 
between g+ and g-. Bottom: resulting density of states (DoS) for a proper choice of the Rabi 
frequencies of the laser beams and their polarization orientation (for details, see [98]). 


non-zero Chern index (we remind that the Chern index for the Lowest Landau Level is 
1). (ii) This lowest band should be very narrow (in units of the recoil energy) to mimic 
the flatness of a Landau level, (iii) It should be separated from the first excited band 
by a large gap, so that the dynamics can safely be restricted to this lowest band even 
in the presence of interactions between particles. We show in fig. 15 the result of an 
optimization for a two-level atom, with an in-plane triangular laser configuration that 
fulfills these three criteria: the lowest band has a Chern index of 1, a width of 0.01 E r 
only, and it is separated from the next band by a gap 40 times larger than its width [98]. 

The notion of an optical flux lattice can be generalized in various ways. First it is 
possible to design optical flux lattice configurations with Chern indices larger than 1. 
Also one can identify situations where: (i) the lowest band has a non-zero Chern number 
(ii) Berry’s curvature (53) calculated for the adiabatic dressed state has a zero flux across 
the unit cell in real space [99]. One obtains in this case a situation that is reminiscent of 
Haldane’s model [65]. Another generalization consists in using an atomic structure with 
more internal states, and designing a time-reversal symmetric configuration correspond¬ 
ing to a Z 2 topological insulator [100]. 
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Time-independent 

Hamiltonian 

Time-dependent 

Hamiltonian (frequency U) 

No use of 

p • A as a Lagrange 

U ~ uj c : rotation 

internal states 

parameter in rotation 

U cj c : lattice shaking 

Using internal 
states 

Berry’s phase 

Spin-orbit coupling 

Laser assisted tunnelling 

Optical flux lattices 

Spin-orbit coupling 


Fig. 16. - Classification of the various procedures used in cold-atom setups to simulate magnetism 
(cyclotron frequency cj c ) or spin-orbit coupling. Schemes in the lower right box have not been 
described in this lecture but are addressed in [101, 102, 80]. 


7. — Conclusion 

We have presented in these notes a series of methods that allow one to simulate, at the 
single-particle level and with neutral atoms, the physics of a charged particle in a static 
external magnetic field. Several of these schemes have been successfully implemented 
in the laboratory and others are currently being investigated. A tentative classification 
of these schemes is presented in fig. 16, where we sort the various procedures along two 
criteria: does it take advantage of the internal atomic structure? Is it based on an 
explicitly time-dependent Hamiltonian? 

It is essential to recall that we described here only a fraction of the activity in this field 
of research. First, as mentioned above, the space limitation did not allow us to discuss 
all the schemes that have been proposed so far at the single-atom level. For example 
we focused here onto the Hamiltonian approach to the problem. Another strategy to 
the general search for ‘topological quantum matter’ is instead to reach these states by 
dissipation: One has to engineer a master equation for which the desired topological 
states are ‘dark’, i.e. they lie in the kernel of this master equation and are reached 
after some relaxation time period (see [103] and refs. in). Second, many other systems 
such as photonics devices are also currently investigated to simulate these magnetic- 
like topological effects. We refer the reader to the recent review articles [104, 105] for 
discussion and references on these photonic implementations of artificial gauge fields. 

We did not address in these notes the role of interactions between particles. These 
interactions are of course crucial if one wants to produce many-body states similar to 
those appearing in the fractional quantum Hall regime. For homogeneous systems, we 
refer the reader to the review [106] where the similarities and differences with respect to 
quantum Hall states are discussed in detail, both in terms of the statistical nature of the 
particles - fermions or bosons and of the type and range of the interaction potential. 
For optical flux lattices, which have properties analogous to the Landau states of the 
free-particle case, one can identify many-body states that are similar to those of bulk 
systems [107]. For optical lattices in the tight-binding regime the situation is much more 
open [108, 109]. Finally let us emphasize that in all our discussions, we considered the 
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gauge fields as static and imposed from the outside. An important challenge for the cold 
atom community is to make these gauge fields dynamical, by coupling their states/values 
to the dynamics of the particles, see e.g. [110, 111] and refs. in. 

* * * 

This work is supported by IFRAF, ANR (ANR-12- 247 BLANAGAFON) and ERC 
(Synergy grant UQUAM). The author thanks Jerome Beugnon, Tom Bienaime, Lauriane 
Chomaz, Nigel Cooper, Laura Corman, Fabrice Gerbier, Nathan Goldman and Sylvain 
Nascimbene for many helpful discussions and remarks on these lecture notes. 


Appendix 1: Landau levels 

Eigenstates with the Landau gauge. - Using the Landau gauge A(r) = Bxu y (cf. 
fig. 17), the Hamiltonian can be written 


(127) 


pr _ Pi . {Py^qBxf 
2 M 2 M 


Since it does not depend on the position operator y but only on p y , we can look for an 
eigenstate basis in the form of plane waves along y, \E^(r) = 'ipk(x) e lky • We suppose 
that the sample has a finite size L y along y and we take periodic boundary conditions 
along this axis. The quantum number k is thus quantized k = (2tt/L y ) n y with n y E Z. 

The function ^k ( x ) 1S a solution of: 


(128) 


h 2 

2 M 


V’fcU) + 


(,hk — qBx) 2 

2 M 


tpk(x) = E ipk(x) 


which can be conveniently written 


(129) 


J^2 -y 

~2M^^ + 2 Mu2c ^ _ Xk ^ ^ k ^ = E ^ k ^ 


with Xk = hk / qB = ki^ . For each fc, this corresponds to the eigenvalue equation for a 
harmonic oscillator centered in Xk, with frequency c j c . We thus recover the Landau level 
spectrum of eq. (12). 

Several remarks are in order at this stage: 

• Although the wave number k along y is generally non-zero, there is no kinetic energy 
h 2 k 2 /2M associated to it. The presence of the magnetic field entails that the only 
influence of this wave number is to shift the center of the fictitious oscillator along 
x. We shall deepen this point later, when we calculate the average velocity for a 
particle in state 

• States belonging to the ground level, also called the lowest Landau level (LLL), 
are obtained by taking linear combinations of the ground-state wavefunctions of 
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Fig. 17. - Left: vector potential for the Landau gauge. Right: probability density of a LLL 
eigenstate of the Hamiltonian, T(r) = r ipk(x)e lky , where 'ipk(x) is a gaussian of width £ mag 
centered in Xk — hk/(qB) = ag . 


the various oscillators. Each ground-sate wavefunction is a gaussian centered in Xk 
with a width f? mag : 

(130) ¥ fc (r) oc e -(*-^) 2 /2<?L g e ife?/ . 


• The various T/Fs are orthogonal because of their different variations with y: 

(131) J n(r) ¥*,(r) d 2 r oc J d y = L v 5 ny , n . y , 

where n y and n' y are the two integers characterizing the quantization of the wave 
numbers k and k'. 


• Two consecutive values of k are separated by 2 7r/L y , which leads to very close 
values for the centers of the corresponding oscillators: 

£ 2 97r 

(132) Xk — = 27r ^ ag if k — k' = —, 

Ly L y 

which is very small compared to ^ ma g if the size L y is large compared to £ mag . 

Finally we can estimate the degeneracy Af of a given Landau level, for example the 
LLL, for a rectangular sample of size L x x L y . First, a given eigenstate will play a 
significant role in the expression of a physical wave packet only if the center Xk of the 
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oscillator is located inside the rectangle. This relates the wave number k and the size 

L'x • 

(133) 0 <x k <L x => 0 <k<-^-. 

mag 

Second, the finite size of L y imposes that k is discrete and a multiple of 2 ^/L y . The 
number of relevant independent states \Et is then: 


(134) 




e 2 

^mag 


27r 


A 

Ott 

z,/l ^mag 


where A = L x L y is the sample area. 

Probability current in a Landau state. - In order to understand why there is no kinetic 
energy term h 2 k 2 /(2 M) in the expression of the total energy of the state \lT, we evaluate 
the average kinetic momentum (II) in this state. Using the expression (13) for II = Mv, 
we find: 

(135) (II) = M(v) = J V* k (r) (-i fiV - gA(r)) V k (r) d 2 r, 

where the expression of 4^ for a state of the LLL is given in eq. (130). 

Clearly the component along x of (II) is zero; indeed A x = 0 and the average of 
the momentum p x in the ground state of the oscillator Hamiltonian is also zero. The 
treatment of the ^/-component is less straightforward. The action of p y on T/.(r) gives a 
factor hk so that 

(136) M(v y (k)) = J $* k (r) (hk - qA y ) V k (r) d 2 r. 

The integration along y can be done explicitly to give 

(137) {v y (k)) = U) c J IV’fcU)! 2 (• x k - x) dx. 

For the ground state of the oscillator centered in Xk, this integral is zero because of 
the symmetry of this state with respect to Xk- Therefore the average velocity and the 
average kinetic momentum are zero in spite of the presence of e lky in the expression of the 
eigenstate: the contribution of the canonical momentum operator is exactly compensated 
by that of the vector potential. 

Eigenstates with the symmetric gauge. - We now choose the symmetric gauge A(r) = 
B(—yu x + xu y )/2 leading to the Hamiltonian 

/\ ^ (jj ^ ^ a ^ 

(138) H = H 0 - y L z ’ with H ° = 2 M + 8 Muj c r2 i L * = %Py ~ VPx- 
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Fig. 18. - Left: Vector potential for the symmetric gauge. Right: probability density for an 
eigenstate ip m in the LLL [cf. eq. (139)], plotted here for m = 12. 


This Hamiltonian is the sum of two operators which commute, so that we can find a 
common eigenbasis to H , H 0 and L z . Let us discuss first separately the two components 
of H. The operator L z is the ^-component of the orbital angular momentum and its 
expression is quite simple in polar coordinates, L z = —i hd^, so that its eigenstates 
are functions of the type F(r) e im( C with m G Z and F an arbitrary function of the 
radial coordinate. The Hamiltonian Ho describes a 2D isotropic harmonic oscillator with 
frequency cj c / 2, and its eigenvalues are (no + l)fwj c /2, no G N. 

To proceed further we have to inject the generic form F(r) e irn(f into the eigenvalue 
equation for Ho and find the functions F(r). The solution is relatively involved in the 
general case (see e.g. [112]), but it becomes quite simple if we restrict to the eigenfunc¬ 
tions of the LLL. Considering the set of functions 

(139) y) = (x + i y) m e~ r / 4 ^ g — r m e irnc f e ~ r / 4 C ag 

for m G N, one can check relatively easily that these functions are both eigenstates 
of Ho with energy (m + l)hoj c /2 and of L z with eigenvalue mh. Therefore they are all 
eigenstates of H with the same energy ftw c / 2, corresponding to the lowest Landau level in 
eq. (12). We plotted in fig. 18 the probability density r|-0 m | 2 . It is rotationally invariant 
around the z axis, it is maximal on a circle of radius r m = a/2 m + 1 £ m ag and its width 
A r ~ f ma g does not depend on m. Note that for a disk shape sample centered at the 
origin, we can recover the LLL degeneracy (134) by counting the number of V’m’s whose 
maximum r m is located inside the disk. 

The general form of a LLL state using the symmetric gauge is a linear combination 
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of the V’m’s with arbitrary amplitudes C m : 

(140) ip(x, y) = C m ip m {x, y) = F(u) e~ r2/4e ™z, 

rri 


where u = x + iy and F(u) = is a polynomial or an analytic function of 

the complex variable u. With this writing, one sees that the restriction to the LLL 
corresponds to an effective passage from the 2D quantum motion in the xy plane to the 
ID motion described by functions of the u variable only. 


Appendix 2: Topology in the square lattice 

In this appendix we intend to give an “elementary” derivation of the expression of 
the Hall conductivity for a square lattice of period a pierced by a uniform flux. More 
precisely we want to show that for an insulator-type filling of some of the subbands of the 
lattice, the conductivity is characterized by an integer number, the Chern index , with the 
following physical meaning: when one applies a force F along one direction of the lattice, 
say x, the Hamiltonian - in the proper frame of reference - is time-periodic with period 
£b = h/(aF). For a zero flux, this is at the origin of the well-known Bloch oscillation 
phenomenon. For a non-zero flux, a current of particles appears along the direction y, 
such that the number of particles crossing a x-oriented link of the lattice (j, l) —>> (j +1, l) 
during the time duration £b is equal to the Chern index. The results derived here can of 
course be found in many instances in the literature (see e.g. [61, 113, 96, 114], and [29] 
for a review), but we thought it might be useful for the reader to have a self-consistent 
derivation of this important result within the present set of lecture notes. 

Band structure and periodicity in reciprocal space. - In this section we consider 
non-interacting particles moving on a square lattice in the tight-binding approximation 
(fig. 19a). In the presence of a uniform flux, we write the single-particle Hamiltonian 
using the Landau gauge A = — Byu x : 

(141) H[a,v] = -jJ2(e~ i2 ^ al+ ^\j + lj)(j,l\ + |j, Z + l)<j, Z|) + h.c. 

3,1 

where we have added the parameter v with respect to the Hamiltonian considered in 
eq. (92). In principle v can be eliminated by redefining the phases of the basis states 
| j, /), but it is convenient to keep it explicitly for the calculations below. 

We restrict ourselves to rational values of <a, a = p'/p with coprime positive 
integers, so that the Hamiltonian is spatially periodic with period a along x and pa along 
y. Let us take the rectangle a x pa as the unit cell in the presence of the magnetic flux. 
This magnetic unit cell is represented in fig. 19b in the case p = 5. A given lattice site 
(j, l) is assigned to the cell labelled by the integers (j, /'), with l = pV + r — 1, r = 1,... ,p. 
We suppose that the sample is formed by N such unit cells, with N integer 1. We 
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Fig. 19. - Relevant regions in real (a-b) and reciprocal (c-d-e) spaces for a square lattice in 
the tight-binding, single mode approximation, (a) Square lattice with a unit cell of size a X a. 
(b) Square lattice in a uniform magnetic field. The flux <f> through a initial unit cell a x a is 
<f> = a0o, where 4>o is the flux quantum. The coefficient a is supposed to be a rational number, 
a — P /P-> so that one recovers the spatial periodicity with an increased unit cell, the so-called 
magnetic cell , with a size a x (pa) for the gauge choice of the text. Here p = 5. (c) First 
Brillouin zone (FBZ) for the square lattice in the absence of magnetic flux, (d) First Brillouin 
zone in the presence of magnetic flux (magnetic Brillouin zone), (e) Reduced magnetic Brillouin 
zone : choosing Bloch vectors in this reduced zone is sufficient to span the whole spectrum of 
the Hamiltonian. 


look for the eigenstates of H[a, v\ in the form of Bloch states: 

(142) \tj}[a,v,q]) = -)= ^ e MoQ*+pl'qy) \u^ v [a, v, q]) 


where the normalized state \uj^ [a, z/, q\) inside the unit cell (j, V) repeats itself in a 
periodic manner over the whole lattice, i.e., it can be written 

(143) I u jy [a,v,q]) =’Y^P r [a,v,q\ e^ r - 1)aq y \j, l = pi' + r - 1), 

r =1 


where the coefficients /3 r are independent of the cell index (j, V). Inserting the expression 
(142) into the equation H\ip) = E |^), we arrive for each set <a, q at the pxp eigenvalue 
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problem for the Hamiltonian in reciprocal space 77 [<a, z/, q\: 


(144) n 1/3 [a, v, q] )) = E[a, v, q] \f3[a, u, q ])). 
Setting 

(145) A r = 2 cos[aq x + 27r((r — l)a + z/)], 
this can be written explicitly 


(146) n = -j 


( Ai e iaq y 0 ... 0 e~ iaq y\ 

e~' aq y A 2 e iaq y ... 0 0 

0 e -i aq y A 3 ... 0 0 

0 0 0 ... A p _i e iaq y 

\ e iaq y 0 0 ... e~ iaq y A p J 


m= 


(Bi\ 

p2 


\pp) 


Here we introduced the notation of curly symbols and double brackets for denoting p- 
component operators and vectors in the reciprocal space (i.e., at a fixed q), in order to 
distinguish them from operators and vectors in the full Hilbert space. 

The eigenvalue system (144-146) corresponds to the problem of a particle moving on 
a ID cyclic chain with p sites, with on site energies —J A r , r = 1,... ,p, and nearest 
neighbor couplings —Je ±iaq y (see fig. 20). It gives rises to p eigenvalues E^[a,j/,q], 
s = 1 setting by convention E W < E^ < ... < E^ p \ For each energy E^ s \ we 

have one eigenvector |/3^)) for the Hamiltonian in reciprocal space hence one Bloch 
vector |-0W) eigenstate of H. 

When q is varied continuously ( 11 ) in the first Brillouin zone (FBZ) associated to the 
magnetic unit cell (the so-called magnetic Brillouin zone , see fig. 19d): 


(147) 


7r 
a 


^ Qx — 

a 


7T 

pa 


<q y < 


7T 

pa ’ 


one finds that these p energies give rise to p (non-overlapping) subbands - called magnetic 
Bloch bands - for each couple [ce,z/]. Actually due to the structure of 77 [see eqs. (145) 
and (146)], it is clear that the E^ : s are not independent functions of q x and z/, but only 
depend on aq x + 2ttu. Hence the band structure obtained when q varies in the FBZ is 
independent of v. 

With the representation in Figure 20, it is clear that the eigenvalue problem for 
77 depends only on the phase exp(±i paq y ). Indeed the only relevant (gauge-invariant) 
parameter for this cyclic chain problem is the total phase accumulated in a round trip 


( n ) For a finite-size sample, the values of q x ,qy are discrete, but they become more and more 
closely spaced as the number of cells N increases. 
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A5 



Fig. 20. - The ID cyclic chain described by the Hamiltonian H of eq. (146) in the case p = 5. 
The on-site energies are —J A r , r — 1,... ,p, and the nearest neighbor couplings are — Je ±iaqy . 


around the chain. Hence this eigenvalue problem is unchanged in the shift q y q y + 
27 t /pa, ensuring that the lower and upper sides of the FBZ of fig. 19d are equivalent. 
Also, since H is invariant in the change q x q x + Ztt/cl, it is clear that the left and right 
sides of the FBZ are equivalent. This equivalence between opposite sides of the FBZ will 
be essential later, when we show that the Chern number must be an integer. 

A remarkable property of the eigenvalue problem of eqs. (144) and (146) is that it is 
also invariant in the change q x —> q x + 27 rp' / (pa) provided one shifts simultaneously all 
indices r —)> r — 1 modulo p (we recall that a = p f /p). This means that a given eigenvalue 
of ft is (at least) p times degenerate, since it can be obtained from p independent Bloch 
vectors [q x ,q y \, [qx + 27T /pa,q y \, [q x + 47T /pa,q y \, .... Therefore if one is interested in 
the spectrum of the Hamiltonian, one can restrict the search by choosing Bloch vectors 
in the reduced magnetic Brillouin zone , where both q x and q y span an interval of width 
27 t /pa (see fig. 19e). 

Constant force and unitary transformation. - We suppose now that a constant force 
F is added along x, resulting for example from a electric field in a Hall-type experiment. 
The single-particle Hamiltonian then reads 

(148) ff t ( ot [a, v, F} = H[a, v] - FX 
where the position operator X in the tight-binding model is 

(149) X = aY J j L700AI- 

3d 


With the form (148) we lose the spatial periodicity of the problem along x. However 
we can recover it, thanks to a unitary transformation generated by the time-dependent 
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operator U(t) = exp(—i XFt/h). The total Hamiltonian after the transformation is 

H tot (t) = U(t) h£> U\t) + ih^MuUt) 

(150) = -jJ2{e- i2 ^ al ^ +t ^\j + l,l)(j,l\ + \j,l + l)(j,l\) + h.c. 

3,1 

where we introduced the Bloch time 

< i5i) * e = s?- 

As announced we recover a spatially periodic problem, with the same unit cell a x (pa) 
as before. The presence of the uniform force F along the x axis is described by the 
additional, time-dependent complex phase 2Tvt/t& for the couplings along the x direction. 

It appears clearly in eq. (150) that the total Hamiltonian is not an independent func¬ 
tion of z/, F and £; the presence of the force F along x is indeed summarized in the 
modification: 

(152) H tot [t,a,vo\ —» H io t[t = 0,a,z/(£)], with v(t) = u 0 + 

tB 

Since the spatial periodicity of the Hamiltonian is conserved, it is still meaningfull to 
consider the Hamiltonian in reciprocal space H. Its structure is identical to eq. (146), 
except that the diagonal coefficients A r are now time-dependent: 

(153) A r = 2 cos {aq x + 27r[(r — 1 )a + i /(£)]} , 

hence a time-periodicity t b- The study of the response of the particles to the force 
F therefore amounts to studying the dynamics associated to the time-periodic p x p 
Hamiltonian H, with (in particular) the possibility for the state of the particle to follow 
adiabatically one of the p eigenstates of this Hamiltonian, when the phase 2ttu entering 
in its diagonal coefficients varies linearly in time [115]. 

Bloch oscillations and adiabatic following. - In the absence of magnetic flux (a = 0), 
the Hamiltonian (150) is separable as a sum of two terms describing the motions along 
x and y. In this case it is well known that the force F causes the phenomenon of Bloch 
oscillations, which takes a particularly simple form in the single-band, tight-binding 
approximation. Suppose that one starts from the Bloch state( 12 ) | if(q)) at time t = 0. In 
the transformed frame, the Hamiltonian H tot (t) is translationally invariant so q remains 
a good quantum number: the state is at any time proportional to | 'ip(q)} and the only 
dynamics is contained in the time-dependent phase acquired by this state. Alternatively, 

( 12 ) For zero flux, the FBZ is (—7r/o, 7t/o] x (—7r/a,7r/a] (fig. 19c) and there is only one Bloch 
state, | fj(q)) — JT i e ia ^ qx+lq F |j, Z), associated to a given quasi momentum q. 



52 


J. Dalibard 


in the initial frame( 13 ) where the Hamiltonian is one finds that the state remains at 
any time a Bloch state with a time-dependent quasi-momentum q(t) moving linearly in 
time across the Brillouin zone at speed q x = F/h; a duration t b is then necessary for the 
particle to travel across the full FBZ and come back to its initial value. First observations 
of Bloch oscillations with atoms in optical lattices were reported in [116, 117] and they 
have found many applications over the last decade, from matter-wave interferometry and 
metrology (see e.g. [118, 119, 120]) to the identification of Dirac points in band structures 
[ 121 , 122 ]. 

The presence of a magnetic flux complicates the matter, but some essential features 
of the Bloch oscillation phenomenon remain valid [71, 72, 73]. Suppose again that at 
the initial time t = 0, the particle is prepared in a state |T(0)) equal to a Bloch state 
\ip( s ^) with quasi-momentum qr, z.e., one of the p eigenstates of H[a , z/o] associated to the 
energy E W [a, z/q, q\. Since the Hamiltonian H tot [a, z/(£)] preserves the spatial periodicity 
a x (pa), the state at time t remains a Bloch state with the same quasi-momentum q. 

In the general case, this Bloch state is a linear combination of the various eigenstates 
of H for this Bloch vector q: 

(154) 7s' (t) \'ip ( ' s ' ) [a,i'{t),q\), v{t) = vo+j~- 

.'=i tB 

If the force F is weak enough (or equivalently £b large enough), the state |T(£)) follows 
quasi-adiabatically the subband s when the parameter v varies in time, so that we will 
have at any time | 7 S | « 1 and \^ s /\ <C 1 for s' ^ s. For the following discussion it will 
be sufficient to give an approximate value for 7 S '(t) at the lowest non-vanishing order. 
We will see that the non-zero Hall current in the y direction originates from the “weak 
contamination” of the state |T(£)) by the Bloch vectors of the subbands s' ^ s. 

An approximate value for 7 ^ can be obtained by a perturbative expansion of the 
time-dependent Schrodinger equation (see e.g. the appendix of [29]): 

(155) 7 s>{t) ~ 7 sit) for s'/s. 

We can express the scalar product in the numerator of eq. (155) in terms of the coefficients 
/3 of the periodic part of the Bloch vector, and take advantage of the fact that they depend 
only on the quantity aq x — 27rz/: 


(156) 


7s'(t) 


i F 7 s (t) 


e( s ') - f;0) 


for s' / s. 


( 13 ) In this case the unitary transformation associated to U is a mere shift of the quasi¬ 
momentum q —>> q + u x Ft/k. 
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The velocity operator and its matrix elements. - Ultimately we want to evaluate the 
average Hall current that appears along the direction y when the force F is applied along 
x. We define the velocity operator V y as the time derivative (in the Heisenberg picture) 
of the position operator Y : 


(157) 


v-^-l 

v dt ~ h 


H,Y 


Y = aY J l\j,l){j,l I- 

3,1 


The expression of V y does not depend on the parameters a, v and F entering in the 
definition of the Hamiltonian and it reads in the tight-binding approximation: 

(158) b = i ^ E (!•?’ 1 + !>0) l\ ~ 1 3, l) O', l + II) • 

3,1 

The velocity operator is invariant in a translation over the square lattice. Hence the 
subspace associated to a given Bloch vector q is globally invariant under the action of 
V y . In other words, for the problem under consideration here, the velocity operator is 
fully characterized by the p x p matrix elements between Bloch states 

(159) (V> (s ' ) [o!,i',<2]| V v \ip {s) [a,v,q]), s, s' = 1,... ,p. 

It is convenient to rewrite this matrix element in terms of the coefficients (3^, [3 <s 1 giving 
the periodic part of the Bloch vector: 

(160) Vy |V’ (S) )=«/3 (S ' ) I Vy |/3 (s) )) 

where we have introduced the p x p matrix V y giving the velocity operator in reciprocal 
space: 

0 -e iaq y 0 ... 0 e- iaq y\ 

e -i aq y Q _ e i aq y _ _ Q Q 

0 e~ iaq y 0 ... 0 0 

0 0 0 ... 0 -e iaq y 

^-e iaq y 0 0 ... e~ [aq y 0 j 

One can immediately check that this velocity operator is connected to the Hamiltonian 
TL by the simple relation 


(161) 


V, .Ja 

Vy = 1 T 


Vy = 


1 &H 
h dq y ’ 


(162) 
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This expression allows us to give an explicit expression for the matrix elements (160) 
that will be useful in the following. Starting from the eigenvalue equation TL\/3^ )) = 
E( s ) |y3( s ) )) 5 taking its derivate with respect to q y and multiplying with ((/ 3^ s ^|, we obtain 

(163) «/3 (s,) | V y |/3 (s) »= l - (E^ -E^)((p^\d qy /3^)) for s ^ s', 

(164) «/3«| V y \^ s) ))^ l -d qy E^- 

The Berry curvature. - Since we now have at our disposal the velocity operator along 
the y direction, we can calculate the average flux along this direction when the force F 
is applied along x and the system is prepared in a given Bloch vector \^ s ^). We suppose 
that the force F is small enough for the adiabatic approximation to hold, so that we can 
use the perturbative expansion of eqs. (154) and (156) for the state of the system |T(£)). 
Using 1 7 S | 2 ~ 1, we obtain 


(165) <*(f)|7|*(f)) 




f E {«^/3 (s) |/3 (s,) ))((/3 (s,) |^/3 (s) ))-c.c.}. 

s'^s 


We can formally add to the sum over s' the contribution of the term s' = s, since this 
term is actually zero( 14 ). Using a closure relation we then obtain 

(166) Vy = (*(f)|7|*(i)) « \d qy E^ + ^ s \q), 

where we have introduced the Berry curvature for a given a and for the subband s 

(167) B^{q) = i ({{ d qx (3^\d qy (3^ )) ({ d qy p {s) \d qx p( s) ))) . 

Eq. (166) is the starting point of the semi-classical study of the dynamics of an electron 
in a magnetic field [29]. It shows that the velocity of the particle includes, in addition to 
the usual group velocity V q E^/H, a second term proportional to the Berry curvature 
B^ of the band [123], which is sometimes called the anomalous velocity [29]. 

From the average velocity V y in eq. (166), we can deduce the average flux through a 
given horizontal link j, / —>> j + 1, l of the lattice. Since a Bloch state corresponds to one 
particle that is delocalized over N lattice cells, with each cell having an area pa 2 , the 
number of particles crossing the considered link of length a in a time interval 5t is 

(168) 


( 14 ) The two quantities (( d qx /3 ^ \/3^ )) and (( d qy /3 ^ |/3^ )) are purely imaginary since |/3^ )) is 
normalized. 
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Conduction from a filled band and Chern number. - The last step in our reasoning is 
to consider a situation where a macroscopic number N) of non-interacting fermionic 
particles are simultaneously present and to evaluate the flux in this case. For simplicity 
we consider a situation with exactly N particles and assume that the gas is initially in its 
ground state, with only the lowest subband 8 = 1 populated and all the other subbands 
empty. We also assume that the subband s = 1 is separated from the next subband 8 = 2 
by a gap( 15 ). Then, for a small enough force F, the population essentially remains in 
the lowest subband at any time. 

Starting from a quantity <F(qr) calculated for one particle in a given Bloch state ^(q), 
we obtain the contribution of the N particles with their Bloch vectors spanning the FBZ 
by the following integral: 


(169) 


*(q) 



dq x dq y 
2tt/ a 2i\jpa 


*(<!)■ 


Here, using the expression of the flux (168) for a single Bloch state, we obtain the flux 
of particles crossing a given horizontal link of length a when the subband 8 = 1 is filled: 


r( s=1 ) 

(170) 5N =- St 

tB 


where we have defined the Chern number associated to the subband s 
(171) C ^ = f[f B^(q)dq x dq y . 

FBZ 

Note that only the second term in the right-hand side expression (166) of V y contributes 
to 5N, since the first term proportional to d qy E^ s \q) has a zero-average over the FBZ. 

The expression (170) leads to the interpretation of the Chern number announced in 
the introduction: when one applies a force F along x, the Bloch oscillation phenomenon 
occurs with the time period te = h/aF. In the presence of a flux through the lattice, 
the Hall current in the y direction is such that £( s=1 ) particles cross a given horizontal 
link of length a during the time interval t b- 

As it is defined, the Chern number clearly depends on the flux a. One could think 
that it is also a function of i/q and F; however C^ is actually independent of these 
quantities. Indeed they enter in the Hamiltonian 7 -i only via the coefficients A r defined 
in eq. (153), through the linear combination aq x — 2i r(z/ 0 — £/£b)- Since one performs an 
integration over q x on the FBZ, the values of i/q and £b (hence F) are irrelevant for the 
value of the integral, hence for the determination of the Chern number. 


( 15 ) This would not hold for a = 1/2 where the two subbands touch at a Dirac point. 



56 


J. Dalibard 


The Chern number is an integer. - In order to show this general property, we first 
introduce Berry’s connection in reciprocal space 

(172) A^(q)=i((p^\V q f3^}), 
which is a vector in the q x , q y plane such that 

(173) B (s) (q) = u z • (V, x A (s) (g)) . 

It is then “tempting” to replace the surface integral (171) of B^ on the FBZ by the 
contour integral of A!> s ^ around the edge of the FBZ: 

(174) = A^ s \q)-dq. 

Z7T Jd FBZ 

However this requires some care since, as emphasized by Kohmoto in [114], eq. (171) 
has a subtle topological nature. First we note that Berry’s curvature (167) is gauge- 
independent, i.e., it is not modified in the change |/3)) —> e l6 \(3 )), where 6 is a smooth 
function of q. On the contrary, Berry’s connection is gauge-dependent so that eq. (174) 
could a priori depend on the gauge choice. Second we notice that if there exists a gauge 
choice which defines a global, single-valued phase of the |/3))’s over the whole FBZ, then 
the contour integral (174) must be zero; indeed the FBZ has the structure of a torus 
(opposite sides correspond to the same physical situation) and its “edge” has thus a zero 
length. 

In the presence of a magnetic field, the situation is made subtle by the fact that it 
is generally not possible to define a smooth global and single-valued phase for the \/3)) ’s 
over the FBZ. When looking for such a definition, a possible strategy could consist in 
setting one of the coefficients of \/3 )) in eq. (146), say the first one /?i, to be real and 
positive everywhere in the FBZ. But this strategy fails if there exists points in the FBZ, 
where this first component vanishes: the phase of |/3)) is ill-defined at these points, which 
introduces singularities in A!> s \ 

Two options for circumventing this difficulty and using eq. (174) have been developed. 
The first one consists in dividing the FBZ in various zones over which the phase of \f3)) 
is separately well-defined, the remaining task being to properly account (via a gauge 
change) for the discontinuities of A!> s ^ at the boundaries between these zones [96]. The 
other option is to relax the condition that \/3)) should be periodic over the Brillouin zone, 
i.e., single-valued over the torus. In this case can be chosen as a smooth function 
over the FBZ [124], and the contour in the integral (174) becomes a ‘true’ rectangle 
(and not a zero-length line on a torus); the challenge in this case is to properly take into 
account the connexion between the choices for the | (3 ))’s on opposite sides of the FBZ. 
This is the strategy that we adopt now. 

The periodicity of the Hamiltonian in reciprocal space entails that its eigenstates on 
two opposite sides of the FBZ differ only by a phase factor. More precisely, there must 
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Fig. 21. - Connecting opposite sides of the FBZ. Due to the periodicity of the problem, the 
states on opposite sides of the FBZ can differ only by a phase factor, which we denote e 1 ^ 1 ^) 
for the upper-lower sides, and for the left-right sides. 


exist a function 4>\{q x ) that relates the lower and upper sides of the FBZ (fig. 21a): 

(175) |/3 W (fe, -)»= eW°W'Hq x , --))) 

pa pa 

and a function (j> 2 (q y ) that relates the left-hand and the right-hand sides of the FBZ 
(fig. 21b): 

(176) / ^ •'/ .„))> 0^1/3^ l,q y )l 

When one evaluates the contour integral (174) following the path ABCD , one can regroup 
the contribution of the segments AB and CD: 


ab J CD 


A {s \q) dq 


A(Qx) d q x = 4 >i(tt/ a) - </>i(-7r/a) 


and the contributions of the segments BC and DA: 

(178) ([ +[ A {s) {q) ■ dq = - [ <j>' 2 (q y ) dq y = 4> 2 (-n/pa) - foin/pa). 

\JbC JDAj J-ir/pa 

The last step in the reasoning is to notice that the states at each corner of the FBZ 
are all equal up to a phase factor. In particular there are two ways to relate the state in 
the lower left corner A and the one in the upper right corner C: 

e W i(*/“) |/3«(_JL))) = e i[<h(«/°)+<h(-*/pa.)] 


a pa 


a pa 


= gi^O/a) \\ = e i[0i(-7r/a)+0 2 (7r/P a )] 1^0) ( _ — ? _ 0L ) \\ 

a 1 pa a 1 pa 


These two ways are equivalent if and only if 


(179) 


/a) + (j) 2 {—^/pa) = 0 i(— 7 T /a) + 02 (tt /pa) modulo 27r. 
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This entails that the sum of the line integrals of A!> s ^ on the four segments AB , BC , CD , 
and DA obtained by adding eqs. (177) and (178) is a multiple of 27r, hence the Chern 
number (174) is an integer. We have proven this result in the specific case of a square 
lattice, but it can be generalized to more complex geometries (see [29] and refs. in). 

One can go one step further and determine the value of the Chern number for a given 
flux a = p'/p and a subband s [61]. We give here the result without proof for the lowest 
subband; the Chern number appears in the solution of the Diophantine equation 

(180) 1 = p'C s=1 +pV 

where V is an integer such that \V\ < p'/2. For the particular case a = 1/p, ie., p' = 1, 
this gives the Chern number C s=1 = 1 (and V = 0). 
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